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FOREWORD 


This  report  was  prepared  by  the  Plax  Corporation  under  USAF 
Contract  No.  AF  33{6l6)-112.  The  contract  was  initiated  under 
Research  and  Development  Order  No.  618-11.  "Improvement  of 
Packaging  Procedures",  and  was  administered  under  the  direction  of 
the  Materials  Laboratory,  Directorate  of  Research,  Wright  Air  De¬ 
velopment  Center,  with  Mr.  M.  E.  Bowman  acting  as  project  engineer. 
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ABSTRACT 


The  temperature  variation  of  the  P -factor  has  been  found  to 
be  accurately  represented  by  the  two  parameter  equation 

P  =  P0e  -Ep/RT 

in  which  Ep  and  PD  are  constants  for  the  particular  permeant 
polymer  system.  These  parameters  are  calculated  from  careful¬ 
ly  constructed  log  vs.  1/T  plots. 

The  size,  shape,  and  polarity  of  the  penetrant  molecule  are 
factors  influencing  Ep  and  for  substances  no  more  polar  than 
acetone  the  expression 

Ep=  0.  0348V  fr  0.  75  V/L*  2.4  44H 

reproduces  the  data  to  within  about  0.  5  kcal/mol.  where  poly¬ 
ethylene  is  the  permeable  film. 

It  is  established  that  for  four  homologous  series  that 

log  PQ=  mEp<4  b 

where  m  aid  b  are  constants  peculiar  to  each  series. 

Master  plots  for  interpolating  homologous  materials  and  extrapolating 
for  temperature  changes  are  included. 

PUBLICATION  REVIEW 

This  report  has  been  reviewed  and  is  approved 
FOR  THE  COMMANDER 


M.  E.  Sorte 
Colonel,  USAF 
Chief,  Materials  Laboratory 
Directorate  of  Research 
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Introduction 


Permeability  has  been  found  by  various  investigations  to  be 
dependent  upon  several  factors  of  which  solubility,  diffusion,  and  temper¬ 
ature  are  important  factors  and  for  a  constant  temperature  is  frequently 
expressed  as 

(1  -  la)  P  -  DS/d  d  =  x  thickness 

(1  -  1  b)  P  -  DS  where  d  r  1 


It  is  an  experimental  fact  that  the  coefficient  of  permeability  as  well 
as  those  of  diffusion  and  solubility  are  often  well  represented  by  the 
following  equations 

-Ep/RT 

(1  -  2a)  P  =  P0e 

-Ed/RT 

(1  -  2b)  D  =  D0e 


(1  -  2c) 


S  -  SQe 


-ES/RT 


to  express  their  dependence  on  temperature.  P,  D  and  S  are  permea¬ 
bility,  diffusion,  and  solubility  expressed  in  grams/24  hour s  /0.  OOlin.  /100 
in.2  of  polyethylene  unless  otherwise  stated.  PQ,  D0,  SQ,  Ep,  Ed,  and 
Es  are  constants  for  a  given  substance.  R  is  the  universal  gas  constant 
equal  to  1.987  calories /degree  Kelvin/mole  and  T  the  absolute  temper¬ 
ature  in  degrees  Kelvin. 


In  this  study  P  ranged  all  the  way  from  immeasurably  small 
values  (less  than  0.  02)  for  several  of  the  hydrogen  bonded  compounds 
at  32°F  to  almost  immeasurably  large  values  for  many  of  the  hydrocarbons 
at  165°F  (more  than  21,000  for  chlorobenzene  at  this  latter  temperature  !): 
altogether,  over  a  million-fold  variation  in  the  P -factor. 


To  check  the  exponential  dependence  of  P  upon  the  inverse  of  the 
absolute  temperature  and  to  evaluate  the  parameters  PQ  and  Ep,  equation 
l-2a  is  rewritten  as 


(1-3)  logP  r  -(Ep/2.  3R)  (1/T)  +  logP0 

showing  that  a  plot  of  the  common  logarithm  of  P  against  the  reciprocal 
of  the  absolute  temperature  should  be  a  straight  line  with  a  slope  of 
-(Ep/2.3R)  and  an  intercept  at  l/T^O  of  logPQ.  This  has  been  done  and 
the  plot  appears  as  Figure  I. 

The  first  step  then  in  the  reduction  of  the  data  is  the  construction 
of  logP  vs.  1/T  plots  as  shown  in  Figure  I.  From  these  the  parameters 
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Fig.  1.  The  Dependence  of  log^F*  upon  1/T. 
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PQ  and  Ep  are  determined  and  compared.  Particular  attention  is  given 
Ep  since  it  is  this  so-called  "activation  energy"  which  determines  the 
temperature  variation  of  the  P -Factor.  Furthermore,  knowing  Ep  it  is 
possible  sometimes  to  estimate  good  values  for  logPQ  and  hence  for  the 
P -Factor  itself. 

Thus  the  accuracy  of  the  logP  vs.  1 /T  plots  is  of  paramount 
importance  in  this  study,  for  the  reliability  of  all  subsequent  calculations 
rests  squarely  upon  the  validity  of  these  plots  and  the  slopes  and  inter¬ 
cepts  obtained  therefrom. 
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2.  The  Master  Data:  Some  Selected  P-Factors  at  32,  70,  100,  130,  and  165°F 


As  a  liquid  diffuses  out  of  a  bottle  initially  filled  to  the  neck, 
the  observed  P-Factor,  which  is  always  a  composite  of  the  liquid  and 
vapor  transmission  rates,  will  generally  decrease  with  time,  for  that 
part  of  the  air  gap  in  contact  with  the  bottle  can  never  be  completely 
saturated  with  vapor  if  the  outside  of  the  bottle  is  "dry".  Hence  the 
vapor-phase  transmission  rate  will  generally  be  somewhat  less  than  the 
liquid-phase  transmission  rate,  and  the  former  becomes  increasingly 
important  as  the  bottle  tends  toward  emptiness,  causing  the  observed, 
instantaneous  P-Factor  to  decrease  with  time.  (Although  when  collapse 
occurs,  this  "tailing  off"  effect  which  must  necessarily  set  in  sometime 
may  be  delayed  until  the  bottle  is  almost  completely  "empty".) 


For  example,  a  bottle  initially  filled  with  100  grams  of  ethyl 
acetate  and  kept  for  nearly  two  weeks  at  165°F  behaved  as  shown  in 
F igur  e  2 . 


Weight  Loss 
(in  grams) 


Fig.  2.  The  Diffusion  of  Ethyl  Acetate  at  165°F 
(Initial  Net  Contents:  100  grams.) 


The  slope  of  the  solid  line  in  Fig.  2  represents  a  P-Factor  of 
1990,  that  of  the  dotted  line  1700.  It  is  the  former,  which  represents 
the  initial  P-Factor  after  steady-state  has  been  established,  that  is  of 
interest  in  this  report.  With  this  in  mind,  the  following  P-Factors  in 
Table  V  -  Part  I  have  been  revised  somewhat  (generally  upwards): 

Ethylene  glycol  monobutyl  ether  at  165°F;  ethyl 
acetate  at  165°F;  amyl  acetate  at  165°F;  allyl 
alcohol  at  70°F,  sec.  -butyl  alcohol  at  100°F  130°F, 
and  165°F;  dibutyl  ether  at  lOO^and  165°F; 
butyraldehyde  at  lOO^Fand  130°F;  aniline  at  165°F; 
nitroethane  at  130°F;  heptane  at  100°F;  tetradecane 
at  70*^  and  100°F;  lubricating  oil  at  130°Fand  165°F; 
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benzene  at  32  °F;  ortho-xvlene  at  70°,  100°,  and 
1 30°F;  para-xylene  at  326,  70°,  100°  and  130°F; 
cyclohexane  at  70°  and  100°F;  chlorobenzene  at  32°, 

70°  and  130°F;  para-chlor otoluene  at  32°,  70°  and 
130°F. 

These  revisions  resulted  in  a  general  improvement  in  the  overall 
appearance  of  the  data  when  the  logarithms  of  the  observed  P-Factors 
were  plotted  against  1/T.  It  should  perhaps  be  noted,  however,  that  such 
revisions  as  have  been  made  do  not  improve  one's  ability  to  predict 
yearly  losses --at  least  not  at  this  point.  In  fact,  as  already  observed, 
in  the  case  of  ethyl  acetate  at  165°F,  the  bottles  are  essentially  empty 
after  14  days.  The  yearly  loss  in  this  case  is  the  same  as  the  14  day 
loss. 

Finally,  with  regard  to  this  reduction  of  the  data,  it  remains  to 
be  mentioned  that  at  32°F  especially,  most  of  the  slow  (and  hygroscopic) 
solvents  exhibited  either  initial  increases  in  weight  or  else  relatively  small 
P-Factors  at  the  start.  Probably  this  is  all  just  a  matter  of  the  degree 
to  which  the  moisture  in  the  air  permeates  inward  at  the  start.  The  effect 
is  most  pronounced  in  the  case  of  the  93%  sulfuric  acid. 
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Table  1 


The  Permeability  of  Polyethylene  to  Forty  Substances  at  32°, 

70°,  100°,  130°  and  L65°F. 

_ Gms/24  hours/0.  001"/100"2  /  Temp. _ 


Substance 

32°F 

f*4 

O 

o 

r- 

100°F 

130°F 

165°F 

Formic  Acid 

0.251 

0.  666 

2.  790 

12.0 

51. 6 

Acetic  Acid 

0.  351 

3.  080 

13.  55 

66.0 

304.  0 

n-Butyric  Acid 

0.  193 

4.  79 

29.  6 

579.  0 

Betnzoic  Acid 

0.  067 

0.  570 

5.  67 

35.  9 

Methyl  Alcohol 

0.253 

1.22 

5.  35 

27.  8 

n-Propyl  Alcohol 

0.072 

0.495 

2.89 

22.4 

168.  0 

Allyl  Alcohol 

0.  161 

0. 65 

4.  06 

23.0 

137.8 

n-Butyl  Alcohol 

0.462 

3.  12 

20.4 

164. 

sec -Butyl  Alcohol 

0.  619 

5.24 

37.9 

280. 

tert.  -Butyl  Alcohol 

0.260 

2.44 

27.  5 

236.  5 

n-Octyl  Alcohol 

0.498 

3.  31 

25.  6 

188.  5 

Phenol 

0.484 

3.28 

24.0 

119.0 

Ethylene  Glycol  Butyl  Ether 

0.  666 

4.40 

34.9 

231. 

Ethyl  Acetate 

1.903 

16.  55 

83.  3 

378. 

1990. 

Amyl  Acetate 

0.  566 

8.  70 

37.7 

269. 

1300. 

Dibutyl  Phthalate 

0.81 

5.7 

27.  75 

Diethyl  Ether 

48.0 

313.0 

Dibutyl  Ether 

407. 

1480. 

4890. 

B  utyr  aldehyde 

0.890 

10.  13 

66.  1 

339. 

Benz  aldehyde 

0.  369 

6.  80 

32.  85 

206.  5 

1060. 

Acetone 

1.  38 

6.  75 

32.45 

184. 

Methyl  Ethyl  Ketone 

3.  67 

12. 6 

60.  0 

326.  5 

1401. 

Diacetone  Alcohol 

0.  308 

1.96 

15.  0 

92.  5 

Acetic  Anhydride 

0.  809 

4.  145 

29.  65 

150. 

Nitroe  thane 

0.  968 

2.  71 

11.4 

61. 3 

290. 

n-Pentane 

96.  9 

526.  0 

iso -Pentane 

47.  8 

270.  0 

Pentene-2 

180. 

695.  0 

n- Heptane 

48.  6 

269.  5 

929. 

2650. 

8160. 

Decane 

9.49 

71.2 

297. 

1220. 

4120. 

Tetradecane 

1.74 

15.  6 

89.  6 

404. 

Benzene 

50.  5 

440. 

1585. 

4480. 

13670. 

Ortho -Xylene 

36.  1 

284.  0 

1167. 

4270. 

16600. 

Par  a- Xylene 

87.  3 

505. 

1760. 

4760. 

16320. 

Cyclohexane 

31. 6 

281. 

1223. 

3730. 

Chlor  obenz  ene 

76.4 

500. 

1630. 

5980. 

21150. 

p  ar  a  -  Chi  or  otolue  ne 

37.  0 

332. 

1357. 

4670. 

18520. 

Dipentene 

15.  0 

128- 

501. 

2030. 

8110. 

Monochloroacetic  Acid 

0.  31 

1.  82 

15. 

59.9 

W  ater 

0.279 

0.  835 

3.94 

18.4 
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3.  The  Master  Plot  and  Its  Uses 


On  a  separate  sheet  the  data  of  Table  1  have  been  plotted  on  semi¬ 
log  paper  as  indicated  by  equation  (1-3).  This  plot  is  of  considerable 
value  in  the  following  respects.  The  Master  Plot  may  be  found  on  page  71  . 

a.  It  may  be  used  to  check  on  the  over -all  consistency  of  the 
data  since  in  a  vast  majority  of  the  cases  the  points  do  fall  very  nearly 
on  perfectly  straight  lines  as  indicated  except  in  those  cases  where 
several  substances  have  nearly  identical  P-Factors.  In  this  respect  one 
observes  that 

i.  On  the  whole,  the  data  show  up  extremely  well;  in¬ 
deed,  very  much  better  than  expected,  for  as  indi¬ 
cated  later,  certain  aspects  of  the  permeation 
process  almost  lead  one  to  expect  departures  from 
linearity.  But 

ii.  Only  aniline  seems  to  show  signs  of .  consistent  curva¬ 
ture.  However, 

iii.  More  often  than  not  the  100°F  point  lies  a  bit  to  the 
right  of  the  best  straight  line  through  the  data.  Also, 

iv.  The  32°F  point  for  methyl  ethyl  ketone  seems  too 
high  by  more  than  a  factor  of  two.  Finally, 

v.  As  indicated  earlier,  the  initial- slope  correction  for 
tailing  off  effects  almost  invariably  brought  otherwise 
seemingly  erratic  data  into  good  straight-line  agreement. 

b.  Thus,  in  view  of  the  well  established  fact  that  the  data  do 
tend  to  fall  almost  exactly  on  straight  lines,  these  plots  may  be  used  to 
"smooth"  data  that  may  happen  to  exhibit  slight  departures  from  linearity. 
Perhaps  more  important  though, 

c.  With  these  plots  one  can  easily  extrapolate  results  obtained 
in  the  neighborhood  of  room  temperature  into  experimentally  difficult-to- 
measure  regions.  For  example,  for  the  alcohols  and  similar  substances 
which  showed  weight  increases  and/or  erratic  behavior  at  32°F,  one  can 
determine  what  the  P-Factors  probably  would  have  been  if  the  relative 
humidity  had  been  zero  both  in  the  32  F"oven"  and  during  weighing  by  ex¬ 
trapolating  down  from  higher  temperatures  where  hygroscopic  effects  were 
neglibible.  Or,  one  can  state  with  confidence  that  the  P -Factor  for 
n-heptane  at  -40°F  and  zero  relative  humidity  is  0.73.  This  is  not  irrmedi- 
ately  obvious  from  the  data  of  Table  2  alone.  Furthermore, 

d.  One  can  accurately  interpolate  with  these  plots.  Thus  the 
P- Factor  for  n-heptane  at  85°F  is  500.  Also, 
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e.  One  can  interpolate  with  respect  to  members  of  a  homologous 
series.  For  example,  the  P -Factor  for  n- octane  at  130°F  is  undoubtably 
within  a  few  percent  of  2000.  With  this  master  plot  one  can  probably  re¬ 
liably  estimate  P -Factors  for  hundreds,  perhaps  even  thousands  of  similar 
compounds.  It  is  important  though  to  classify  a  compound  according  to  its 
dominant  functional  group.  Thus  diacetone  alcohol,  CHjCOCH^COHfCH^^ 
belongs  with  the  alcohols  rather  than  the  ketones  as  would  be  the  case  if 
the  hydroxyl  group,  which  takes  precedence  over  all  others  studied,  were 
absent. 

Likewise,  it  should  be  noted  that  chlorine  has  no  drastic  effect  upon 
P-Factors,  behaving  in  this  respect  very  much  like  a  methyl  group.  This 
is  somewhat  of  an  oversimplification,  but  observe  that  chlorobenzene  has 
permeability  characteristics  similar  to  benzene  itself,  paxa-chlorotoluene  to 
par  a -xylene,  and  chloroacetic  acid  to  acetic  acid. 

As  a  further  example  of  the  importance  of  properly  assessing  the  ef¬ 
fects  of  functional  groups,  the  following  (with  editorial  remarks  in  parenthe¬ 
ses)  is  quoted  from  an  article  "Polyethylene  Packaging  Problems"  by  Wight, 
Tomlinson  and  Kirmeier,  Drug  and  Cosmetic  Industry,  72,  767  (1953). 

"Resinous  and  crystalline  materials  familiar  to  the  trade 
as  heavy,  persistent,  and  of  low  volatility  are  not  lost  in  ap¬ 
preciable  amounts.  Alcohols  are  superior  to  the  correspond¬ 
ing  acetates  by  a  factor  ranging  from  two  to  ten,  and,  likewise, 
but  to  a  smaller  degree,  are  less  permeable  than  the  respective 
aldehydes;  (Superior  means  smaller  P-Factor.)  dimethyl  acetals 
are  unexpectedly  poorer  than  the  parent  aldehydes  (that  is, 
they  have  larger  P-Factors).  Terpenes  and  terpenic  oils,  low 
hydro-carbons  and  low  ketones  and  esters  are  phenomenally  bad 
while  hydroxycitronellal  (we  repeat,  hydroxycitronellal)  is  practi¬ 
cally  non-permeable.  Others  skip  around  in  a  most  unpredicta¬ 
ble  maimer,  eugenol  diffusing  quite  slowly,  anethole  very  rapidly.  " 

Insofar  as  these  observations  overlap  those  of  the  present  study,  the 
order  is  the  same;  namely  (in  order  of  increasing  P-Factor): 

Alcohols 

Acids 

Nitr  oderivative  s 

Aldehydes 

Ketones 

Esters 

Ethers 

Hydrocarbons 

(This  is  shown  in  Table  2.) 
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Table  2 


ORDER  OF  P -FACTOR  IN  RELATION  TO  FUNCTIONAL  GROUP 
{From  Table  V  Part  I  of  This  Report) 


Group 

Specific  Material 

P- Factor 
gms/24  hr. 
0. 001"/100"2 

Alcohols  -C-OH 

Methyl  alcohol 

1.22 

n  Propyl  alcohol 

0.49 

Acids  -  C»0  O  H 

Acetic  acid 

3.08 

Butyric  acid 

4.  79 

Nitro  derivatives 

Nitroethane 

2.71 

-C-N02 

Nitrobenzene 

4.93 

Aldehyde  s  -  He  *0 

Acetaldehyde 

6.05 

Ketones  C  -  O 

Acetone 

6.  75 

Methyl  ethyl  ketone 

12. 60 

Esters  - C'^-O-R 

Ethyl  acetate 

16.  55 

Amyl  acetate 

8.75 

Ethers  -C-O-C 

Ethyl  ether 

313. 

Butyl  ether 

85. 

Hydrocarbons 

n-Pentane 

526. 

Benzene 

440. 
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The  reason  why  acetals  permeate  more  rapidly  than  the  correspond¬ 
ing  aldehydes  is  that  the  exposed  oxygen  of  an  aldehyde  becomes  an  ether 
linkage  in  the  corresponding  acetal,  and,  as  noted  above,  ethers  permeate 
more  rapidly  than  aldehydes.  The  increase,  however,  will  not  be  as  great 
as  one  might  expect  from  this  factor  alone  since  the  acetals  are  bulkier 
than  the  original  aldehydes.  (The  effect  of  bulkiness  upon  P-Factors  will 
be  demonstrated  when  we  compare  the  slopes  of  the  curves  LogP  vs.  1 /T 
for  ter -butyl  alcohol  and  n-butyl  alcohol,  ortho-xylene  and  para-xylene,  and 
cyclohexane  and  pentane  or  heptane.) 

As  for  hydroxycitronellal,  anytime  you  have  a  hydroxyl  group  present 
you  can  to  a  first  approximation  just  about  forget  about  the  rest  of  the 
molecule.  The  permeability  through  polyethylene  will  be  low.  This  then 
explains  eugenol  and  anethole.  The  former  is  essentially  an  alcohol,  the 
latter  an  ether  (Figure  3-1). 


H 

0 

OC.Li; 

/\  ' 

\/ 

\  / 

CU_0!  =  t*._ 

Cr!  -  CKCh  , 

Eugenol  Anethole 

Fig.  3.  The  Structures  of  Eugenol  and  Anethole. 


4.  Determination  of  the  Parameters  E  and  P  . 
_ P _ o_ 

Ep  values  were  calculated  from  the  slopes  of  curves  similar  to  those 
presented  in  the  Master  Plot  using  either  the  1/T=2.80  and  l/T-3.85  inter¬ 
cepts  (Method  I)  or,  for  the  slower  permeators,  the  1/T;2.80  and  1/T=3.40 
intercepts  (Method  II  ): 

(4-la)  E^  —  4,360  [logP(  1 /Tt 2 . 80)  -  logP(  1 /T=3.  85)j  ,  or 

(4-lb)  Ep  ~  7,  626  [logP(l/T=2. 80)  -  logP(l /T=3. 40)] 

calories  per  mole. 

In  all  cases  PD  was  calculated  from 

(4-2)  logP0=  logP(l/T»2.80)  +  0.6ll9(Ep/1000)  - 

using  the  Ep  value  obtained  from  either  method  I  or  II  above. 
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The  intercept  data  and  results  of  the  calculations  axe  tabulated  in 
Table  3. 

In  every  case  refers  to  the  1/T=2.80  intercept  and  Pj  to  the 
1/T=3.85  or  1/T=3.40  intercept  depending  upon  which  method- -equation 
(4- la)  or  (4- lb)  respectively- -was  used.  This  is  indicated  in  the  second 
column. 

The  Ep  values  were  calculated  to  five  significant  figures  and  then 
rounded  off  to  four  although  they  are  at  best  reliable  to  only  three  signifi¬ 
cant  figures.  The  PQ  values  have  been  rounded  off  to  three  significant 
figures,  again  undoubtably  one  more  than  is  actually  warranted.  In  the 
case  of  the  pentanes  the  calculations  are  based  upon  only  two  points. 
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Table  3 


The  Parameters  Ep  and  PQ  in  the  Equation  PxP0e“^p^^ 
(Listed  in  order  of  decreasing  Ep) 


Substance 

Method  P2* 

Pi  Ep 

kc  al/mol 

logPo 

0 

*— » 

O 

1 

O) 

ter -Butyl  Alcohol 

n 

710 

0.  260 

26.  21 

18. 886 

7,  690, 000 

sec -Butyl  Alcohol 

II 

650 

0.  619 

23.04 

16.  911 

81, 500 

Butyl  Carbitol 

II 

620 

0.  680 

22.  57 

16. 604 

40,200 

n- Octyl  Alcohol 

II 

445 

0.  500 

22.49 

16.411 

25,  800 

n- Butyl  Alcohol 

II 

395 

0.462 

22.  36 

16.278 

19,000 

iso-Propyl  Amine 

I 

33, 900 

0.  260 

22.  30 

18. 177 

1, 500, 000 

n-Propyl  Alcohol 

II 

381 

0.48 

22.  11 

16. 112 

12, 900 

Di acetone  Alcohol 

n 

357 

0.  308 

21. 91 

15.769 

5,  870 

Benzoic  Acid 

11 

83 

0.  12 

21. 66 

15.  171 

1,480 

Phenol 

11 

293 

0.484 

21.22 

15.449 

2,810 

Acetic  Anhydride 

11 

357 

0.  809 

20.  17 

14. 894 

783 

Allyl  Alcohol 

n 

304 

0.  680 

20.21 

14. 851 

709 

Dibutyl  Phthalate 

n 

66 

0.  153 

20.  09 

14. 114 

130 

Monochloroacetic  Acid  II 

155 

0.  39 

19.82 

14. 318 

208 

Amyl  Acetate 

1 

2,820 

0.  100 

19.40 

15.  323 

2,  100 

Benz  aldehyde 

11 

2,  360 

6.  80 

19.  37 

15.227 

1,690 

B  utyr  aldehyde 

1 

3,400 

0.  145 

19.05 

15.  191 

1, 550 

Acetone 

n 

1,700 

6.  30 

18.  54 

14.724 

530 

Water 

H 

39.4 

0.  150 

18.45 

12.886 

7.  69 

n -Butyric  Acid 

11 

1,290 

5.  10 

18.  32 

14. 324 

211 

Methyl  Ethyl  Ketone 

1 

2, 830 

0.207 

18.  03 

14.486 

306 

Tetradecane 

1 

3,  800 

0.284 

17.99 

14. 588 

388 

Nitroe  thane 

n 

357 

0.  809 

17.  88 

13.721 

52.6 

Acetic  Acid 

11 

644 

3.  00 

17.  78 

13.  690 

48.  9 

Methyl  Alcohol 

11 

241 

1.  16 

17.  67 

13. 197 

15.7 

Ethyl  Acetate 

11 

3,  670 

17.4 

17.72 

14.408 

256 

Formic  Acid 

n 

100 

0.  65 

16.  68 

12.206 

1.  61 

Cyclohexane 

1 

32,200 

6.  3 

16.  17 

14.402 

252 

Dipentene 

1 

15,  700 

3.  35 

16.  00 

13. 989 

95.  7 

Dibutyl  Ether 

1 

10, 900 

2.41 

15.  94 

13.789 

61.6 

p  -  Chlor  otoluene 

1 

35, 000 

8.  15 

15.  84 

14.236 

172 

0- Xylene 

1 

30, 300 

8.  00 

15.  60 

14.  028 

106 

n- Decane 

1 

7,  650 

2.06 

15.  56 

13.407 

25.  6 

Chlorobenzene 

1 

35, 500 

19.8 

14.  18 

13.230 

17.0 

p- Xylene 

1 

28,  300 

24.  3 

13.  37 

12. 632 

4.29 

Benzene 

1 

23,  800 

21.  5 

13.27 

12.498 

3.  15 

n- Heptane 

1 

14, 000 

13.  0 

13.22 

12.235 

1.72 

iso-Pentane 

1 

14, 000 

13.4 

13.  16 

12.200 

1. 58 

n-Pentane 

1 

25, 000 

29.  0 

12.80 

12.230 

1. 70 

*PD,  Pj  and  P2  have  the  same  units  as  P-Factor  -  gms/24  hrs/0.  001 "  /100  in2 
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The  Physical  Significance  of  the  Activation  Energy  for  Permeation,  E  . 


It  has  been  argued  in  considerable  detail  that  insofar  as  diffusion  in 
polyethylene  is  Fickian, 


(5-1) 


Ep~  Ed  +  Es 


where  Es  is  essentially  the  heat  of  solution  of  the  permeant  in  polyethylene 
and  E<j  is  the  energy  of  activation  for  diffusion.  . 


It  has  been  established  that 


(5-2)  P  =  DS 

where  D  is  the  diffusion  constant.  From  the  comparison  of  propanol  and 
decane  it  seems  probable  that  Eg  is  proportional  to 

(5-3)  AAH  =•  AH  .  -&H  .  ..  ,  .  , 

vaporization  vaporization  of  a  hydro- 

of  permeant  carbon  of  similar  size. 

This  quantity &&H  will  be  zero  for  hydrocarbons  such  as  decane  which 
are  non-polar.  Therefore  Eg  will  be  small,  S  large,  and  the  P -Factor, 

P  large  in  such  cases. 

For  polar  compounds,  however,  we  anticipate  that  will  be  in¬ 

herently  positive.  Thus,  using  the  group  contributions  to  the  cohesive 
energy  of  organic  liquids  as  given  in  Table  4,  we  calculate  that  for  propanol 

aah  =  -ah^3ch2cH3  =(i.s  + 2*1.0  +  7.3)  -  (2*1.8+  1.0) 

or  6.  5  kcal/mole.  And  this  presumably  is  just  about  an  upper  limit  to  Eg 
(Indeed,  the  absorption  tests'  data  to  be  discussed  later  suggest  that  the 
heat  of  solution  of  propanol  in  polyethylene  is  around  4.  7-4.  8  kcal/mole.) 
Since  Ep  itself  is  22.  1  kcal/mole,  this  means  that  for  propanol  E^  is  at 
least  22.  1  -  6.  5  =  15.  6  which  is  several  times  greater  than  Eg.  Further¬ 
more,  insofar  as  Eg  s  O  for  the  hydrocarbons,  E^  =  Ep  which  is  greater 
than  12  kcal/mole  in  all  cases  investigated  in  this  project. 
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Table  4 


GROUP  CONTRIBUTIONS  TO  THE  COHESION  ENERGY 


OF  ORGANIC 

LIQUIDS  * 

Group 

Cohesion  in 
kcal  /mole 

Group 

Cohesion  in 
kcal /mole 

-ch2- 

1.0 

-CO 

4.  3 

=  CH- 

1.0 

-CHO 

4.7 

-0- 

1.  6 

-cooch3 

5.  6 

rCH2 

1.8 

-cooc2h5 

6.2 

-ch3 

1.8 

-no2 

7.2 

-F 

2.  1 

-OH 

7.  3 

-Cl 

3.4 

-COOH 

9.0 

-Br 

4.3 

-conh2 

13.2 

-I 

5.  0 

-CONH- 

16.2 

*  Data  from  Dunkel  (Z.  physik.  Chem.  A  138,42  (1928) 


In  fact,  the  argument  that  led  to  the  introduction  of  quantity  ^>*-H 
compels  one  to  ascribe  such  variations  as  there  are  in  the  Ep  values  for 
the  hydrocarbons  primarily  to  E^.  And  these  variations  are  not  insignifi¬ 
cant,  Ep  ranging  all  the  way  from  12.8  for  n-pentane  to  18.0  for  tetradecane. 

This  all  seems  to  say  that  D  rather  than  S  is  the  governing  factor  in 
the  product  DS.  In  a  sense  this  is  correct.  Our  earlier  conclusions  re¬ 
garding  the  dependence  of  P  upon  S  and  hence  upon^^H  are,  however, 
still  valid  provided  we  supplement  them  with  the  following  observations. 

a.  As  one  might  expect,  size  and  shape  play  an  important  role 
in  the  permeation  process.  They  affect  Ep  primarily  through  their  influ¬ 
ence  upon  E^,  the  kinetic  part  of  Ep,  having  relatively  little  to  do  with  Eg, 
the  static  or  equilibrium  part  of  E  .  Thus  whereas  the  solubilities  of  the 
two  xylenes,  ortho  and  para,  are  Nearly  identical,  their  activation  energies 
for  permeation  are  significantly  different:  15.6-  and  13.4  respectively. 

This  difference  we  ascribe  to  the  bulkiness  of  the  ortho  molecule,  for  both 
isomers  have  essentially  the  same  molar  volume. 
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Or  take  the  case  of  ter-  and  n-butyl  alcohol,  again  rather  two  differ¬ 
ently  shaped  isomers  with  nearly  the  same  molar  volumes.  Their  Ep  values 
are  26.2  and  22.4  respectively,  yet  the  absorption  tests  indicate  that  their 
heats  of  solution  differ  by  less  than  0.  1  kcal/mole.  The  remaining  differ¬ 
ence  must  be  due  to  the  bulkiness  of  the  tertiary  construction  as  seen  during 
diffusion.  Finally,  when  we  consider  the  normal  aliphatic  hydrocarbons,  we 
find  that  their  solubilities  are,  relatively  speaking,  all  just  about  the  same, 
yet  there  is  a  steady  increase  in  Ep  in  going  from  pentane  through  heptane 
and  decane  to  tetradecane.  In  this  case  the  molecules  are  all  the  same 
general  shape  but  of  rather  different  sizes,  and  in  general  the  larger  the 
molecule  the  more  difficult  its  motion  during  diffusion  from  one  equilibrium 
position  to  the  next.  Here  again  the  effect  upon  Ep  is  largely  through  E,j 
rather  than  Eg. 

b.  Nevertheless,  other  things  being  the  same --namely  size  and 
shape- -polarity  as  measured  by&AH  (which  depends  upon  AH^  rization^ is 
the  determining  and  distinguishing  factor,  and  in  the  case  of  tne  more  polar 
compounds  it  easily  dominates  factors  of  size  and  shape.  Thus  the  addition 
of  CH.2  groups  to  pentane  increases  Ep  by  only  5.2  kcal/mole  whereas  the 
replacement  of  one  hydrogen  atom  by  a  single  OH  group  increases  Ep  by 
more  than  8  kcal/mole  ! 

c.  However,  and  this  is  where  we  modify  our  earlier  conclusions 
somewhat,  the  full  effect  of  polarity  upon  Ep  may  not  be  felt  at  the  Eg 
stage  owing  to 

i.  an  appreciable  interaction  of  polyethylene  with  polar 
sub  stanc  e  s  and 

ii.  solvent  association  effects  in  the  polymer  systems. 

Evidence  supporting  these  conclusions  will  be  presented 
in  later  sections  when  we  examine  in  detail  the  structure 
of  Ep,  the  relation  between  Ep  and  PQ  ,  and  the  temper¬ 
ature  independent  part  of  the  polymer  interaction 
parameter  as  revealed  by  the  absorption  tests. 

d.  Nevertheless,  that  part  of  &AH  which  fails  to  appear  in  Eg  owing 
to  polymer  interaction  and  solvent  association  shows  up  eventually  in  Ep 
through  the  E^  term  except  in  extreme  cases  of  polarity  as  in  the  hydrogen 
bonded  alcohols  and  acids  where  it  now  seems  probable  that  a  dimer  rather 
than  the  monomer  is  the  diffusing  entity. 

e.  To  summarize  then, 

i.  Ep  (and  therefore  P)  is  very  structure  sensitive,  de¬ 
pending  upon  the  size,  shape,  and  polarity  of  the  permeant 
molecule. 

ii.  But  whereas  Eg  and  Ep  both  increase  with  increasing 
polarity  as  measured  by  J*AH, 
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iii.  Only  E,  is  markedly  dependent  upon  the  size  and 
shape  factors. 

iv.  However,  when  it  exists,  the  polarity  factor  may 
greatly  outweigh  all  others  in  its  influence  upon  E^. 

v.  Therefore  S  may  be  used  as  a  qualitative  guide  to 
P,  although  it  is  perhaps  important  to  note  that 

vi.  The  variations  in  S  withAAH  will  not  be  as  great 
as  those  in  P  itself  since  D  also  varies  withAdH. 

Also, 

vii.  Eg  (and  therefore  S)  will  not  properly  reflect  the 

effect  of  variations  in  size  and  shape  upon  E^  (or  P). 

So  much  for  the  general  structure  of  E  .  We  turn  next  to  a  consider¬ 
ation  of  special  cases. 

Now  it  would  be  nice  if  we  could  by  a  consideration  of  special  cases  es¬ 
tablish  in  a  general  way  the  dependence  of  both  Eg  and  E^  individually  upon 
each  of  the  factors:  size,  shape,  and  polarity.  Then  by  synthesizing  the  two, 
we  would  have  an  elegant  derivation  of  a  general  expression  for  Ep.  This 
procedure  presupposes  that  we  have  some  method  for  determining  either  Eg  or 
E^  in  each  case,  and  it  was  hoped  initially  that  perhaps  the  temperature  vari¬ 
ation  of  the  solubility  would  permit  a  determination  of  one  of  these,  namely  Eg. 
E<j  could  then  be  obtained  by  difference  and  the  two  studied  separately.  Un¬ 
fortunately  the  solubility  behavior  in  a  number  of  crucial  instances  appears 
anomolous  from  the  point  of  view  of  present  polymeric  theory.  Probably  the 
partial  crystallinity  of  polyethylene  is  responsible  for  the  observed  departures 
from  theory.  In  any  event,  a  reliable  decomposition  of  E  into  two  terms,  one 
associated  with  solubility  and  the  other  with  diffusion,  has^probably  not  been 
effected.  Experimentally  there  still  remains  the  possibility  of  independently 
determining  E^  from  time  dependent  absorption  and  desorption  studies.  How¬ 
ever,  lacking  such  information,  we  turn  now  to  E^  itself. 

6.  The  Qualitative  Dependence  of  E^  Upon  Size,  Shape,  and  Polarity. 

Water  and  tetradecane  form  an  interesting  comparison.  Molecule  for 
molecule,  tetradecane  is  more  than  140  times  larger  than  water,  yet  the 
activation  energy  for  permeation  of  the  hydrocarbon,  CH^C^C^C^CH^ 

CH2CH2CH2CH2CH2CH2CH2CH2CH3,  is  18.0  kcal/mole  whereas  that  for  HzO 
is  18.4  kcal/mole! 


Comparisons  such  as  this  are  presented  in  a  more  systematic  manner 
in  Table  5. 
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Table  5 


Ep  and  Size,  Shape,  and  Polarity. 


Substance 

Molar  Volume 

Ad  H 
kcal/mol 

EP 

kcal/mol 

Pentane 

115 

0.0 

12.8 

Decane 

194 

0.0 

15.  6 

Tetradecane 

259 

0.0 

18.  0 

Ethyl  Acetate 

98 

2.9 

17.  7 

Amyl  Acetate 

148 

2.  5 

19.4 

p- Xylene 

124 

0.4 

13.4 

o- Xylene 

121 

0.  6 

15.6 

n-Pentane 

115 

0.  0 

12.  8 

Cyclohexane 

108 

0.  5 

16.  2 

n-Butyl  Alcohol 

94 

5.4 

22.4 

ter -Butyl  Alcohol 

92 

4.8 

26.2 

Benzene 

89 

0.  6 

12.  5 

Butyr  aldehyde 

88 

2.4 

15.2 

Phenol 

88 

5.2 

21.2 

Heptane 

146 

0.  0 

13.2 

Amyl  Acetate 

148 

2.5 

19.4 

Acetone 

74 

3.  0 

18.  5 

n-Propyl  Alcohol 

75 

6.  3 

22.  1 

In  the  pentane -decane -tetradecane  series  the  shape  and  polarity  factors 
are  constant  so  that  one  is  observing  here  the  effect  upon  Ep  of  size  alone. 
This  is  also  nearly  true  of  the  two  acetates. 

• 

On  the  other  hand  the  shape  factor  differs  for  the  two  xylenes.  Their 
size  and  polarity  factors  are  nearly  identical.  Likewise,  the  principle 
difference  between  n-pentane  and  cyclohexane  or  n-butyl  and  ter -butyl  alco¬ 
hol  is  one  of  shape. 

In  the  remaining  examples  of  Table  5  the  polarity  factor  varies,  the 
size  and  shape  within  each  set  being  essentially  constant  (with  the  exception 
of  butyr aldehyde).  Observe  that  unless  the  molecules  differ  widely  in 
geometry,  the  previously  mentioned  order  of  decreasing  Ep  (alcohol, 
aldehyde,  ketone,  ester,  and  hydrocarbon)  seems  to  be  preserved. 
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The  problem  before  us  then,  is  one  of  reducing  these  features  of 
Table  5  to  some  sort  of  analytical  statement  that  could  be  used  to  calculate 
reasonably  reliable  Ep  values  for  substances  of  unknown  permeability  from 
readily  available  data.  Such  data  for  the  substances  under  consideration  in 
this  report  are  assembled  in  Table  6. 


Column  three  gives  the  vapor  pressures  at  25°C  estimated  from 
Dreisbach's  "P-V-T  Relationships  of  Organic  Compounds"  which  has  been 
chosen  as  a  convenient  source  of  thermochemical  data.  Column  four  lists 
approximate  heats  of  vaporization,  uncorrected  for  non-ideal  gas  behavior, 
as  estimated  from  Dreisbach.  Then  follows  aaH  defined  in  equation  (5-3) 
and  evaluated  from  Dreisbach's  heats  of  vaporization;  aaH  evaluated  using 
Table  4  (listed  under  AAH);  the  molar  volume  in  cc  at  room  temperature; 
and  L,  the  length  in  centimeters  of  the  Taylor -Hirschfelder- Fisher 
molecular  models. 
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Table  6 


Supplementary  Data  Sheet:  Size.  Shape,  and  Polarity  Factors 


Substance 

B.P. 

V.P. 

AH  _ 

AAH 

AAH’ Volume 

L 

Vol/L 

•c 

mm  /Hg 

k/car/mol 

cms^/mol 

cms* 

ter -Butyl  Alcohol 

82.  8 

45. 

10.  5 

4.  8 

7.  3 

94.  0 

5.9 

15.  9 

sec -Butyl  Alcohol 

100. 

20. 

11.  1 

5.4 

7.  3 

91.  7 

7.  7 

11.9 

n- Octyl  Alcohol 

195. 

0.2 

14.  3 

4.8 

7.  3 

158.  0 

14.  6 

10.8 

n-Butyl  Alcohol 

117.  7 

6.  4 

11.7 

6.  0 

7.  3 

91.  6 

9.  1 

10.  0 

iso -Propyl  Amine 

34.0 

550. 

6.  7 

2.  1 

— 

91.0 

6.7 

13.  6 

n-Propyl  Alcohol 

97.  2 

20.  1 

11.0 

6.  3 

7.  3 

74.  8 

7.  8 

9-  6 

Diacetone  Alcohol 

165. 

1.0 

13.  3 

4.  7 

— 

125.  0 

9.  1 

13.  8 

Phenol 

182. 

0.7 

12. 7 

5.2 

7.  3 

87.  7 

8.  5 

10.  3 

Allyl  Alcohol 

96. 

25. 

10.  9 

6.  1 

7.  3 

67.8 

7.  7 

8.  8 

Chloroacetic  Acid 

189. 

0.  3 

14.  1 

9.4 

10.  6 

59.7 

7.  3 

8.2 

Amyl  Acetate 

148.* 

5. 

.10.  9 

2.  5 

2.4 

148.  0 

13.2 

11.  2 

Benz  aldehyde 

179.  5 

1.  5 

10.  9 

2.  5 

— 

101. 

9.0 

11.2 

B  utyr  aldehyde 

75.7 

100. 

8.  1 

2.4 

3.  6 

88.  2 

8.  5 

10.4 

Acetone 

56.  5 

230. 

7.7 

3.0 

3.  3 

73.  8 

6.  5 

11.  3 

Water 

100. 

23.  5 

10.  5 

— 

— 

18.  0 

3.7 

5.2 

nButyric  Acid 

164. 

1. 

13.  1 

7.4 

9.0 

92.  0 

10.  6 

8.  7 

Methyl  Ethyl  Ketone 

79.  6 

105. 

8.  5 

2.  8 

3.  3 

89-  5 

7.9 

11.  3 

Tetradecane 

252.  5 

0.  0 

14.4 

0.  0 

0.  0 

259-  0 

21.  6 

12.  0 

Nitroethane 

114.  8 

20. 

9.  5 

6. 

6. 

71.  3 

6.0 

11.9 

Acetic  Acid 

118.  1 

17. 

10.0 

6.  1 

9.0 

57.  3 

6.  3 

9.  1 

Methyl  Alcohol 

64.  6 

100. 

9.9 

7.  1 

7.  3 

39.  5 

4.  7 

9.  1 

Ethyl  Acetate 

77.2 

100. 

8.  5 

2.  8 

2.4 

97.9 

9.2  . 

10.  6 

Formic  Acid 

100.  7 

46. 

8.  7 

6. 

9. 

37.  5 

5.  6 

6.  7 

Cyclohexane 

81.4 

100. 

8.  1 

0.  6 

0. 

81.4 

6.  9 

15.  6 

Dipentene 

177. 

2. 

11.  5 

0.  0 

0. 

162. 

10.  7 

15.  1 

Dibutyl  Ether 

142. 

7. 

10.  3 

0.  8 

0.  6 

169.  5 

15.  1 

11.  3 

p  -  Chlor  otoluene 

162. 

3. 

10.  8 

1.  6 

1.  6 

118. 

9.  5 

12.4 

o- Xylene 

144. 

17. 

10.  1 

0.  6 

0. 

121. 

8.7 

13.9 

p- Xylene 

138.  5 

26. 

9.9 

0.4 

— 

— 

9.7 

12.  7 

n-Decane 

174. 

1.  5 

11.4 

0.  0 

0.  0 

194. 

16.2 

12.  0 

Chlorobenzene 

132. 

12.  5 

9.  8 

1.4 

1.  6 

101.  5 

8.  6 

11.  8 

Benzene 

80.  0 

100. 

8.  1 

0.  6 

0. 

88.  9 

7.  0 

11.7 

n- Heptane 

95.  5 

50. 

8.  4 

0.  0 

0.  0 

146.  5 

12.2 

12.  0 

n-Pentane 

36.2 

500. 

6.  6 

0.  0 

0.0 

115.2 

9.  5 

12.  1 

In  constructing  this  table,  the  following  values  have  been  used  for  the 
heats  of  vaporization  of  the  normal  aliphatic  hydrocarbons:  ethane,  approx. 
3.8;  propane,  approx.  4.7;  butane,  5.7;  pentane,  6.6;  hexane,  7.6;  heptane, 
8.6;  octane,  9.6;  nonane,  10.6;  decane,  11.6;  undecane,  12.6;  dodecane, 
13.6;  tridecane,  14.6  except  that  tetradecane  was  taken  as  determined  from 
Dreisbach. 


*  at  737  millimeters  of  mercury. 
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The  substances  in  Table  6  have  been  listed  as  in  Table  3  according  to 
decreasing  Ep.  Observe  that  the  progressive  diminuition  of  Ep  does  not  seem 
to  be  reflected  by  any  pronounced  trend  in  either  the  boiling  points  (given  in 
degrees  centigrade),  the  vapor  pressures  (approximate  values  in  mm),  the 
heats  of  vaporization  (approximate  values  in  kcal/mole),  the  molar  volumes 
(expressed  in  cm^).  the  lengths  L  of  the  Fisher-Taylor-Hirschfelder  molecu¬ 
lar  models  (measured  in  cm)  or  the  ratios  of  the  molar  volumes  to  L.  '  If 
there  is  any  significant  correlation  at  all  it  is  with^H,  illustrating  as 
mentioned  many  times  previously  that  polarity  plays  a  prominent  role  in 
permeability  behavior.  It  is  for  this  reason  that  equilibrium  solubility  studies, 
which  reflect  polarity  more  than  any  other  single  factor,  furnish  considerable 
information  about  P- Factors. 

In  the  next  section  we  begin  a  quantitative  study  of  the  dependence  of 
Ep  upon  size,  shape,  and  polarity. 

7.  The  Dependence  of  Ep  Upon  Size. 

The  strategy  here  is  essentially  that  indicated  by  Table  5:  to  study  the 
variations  in  E_  for  a  series  of  compounds  differing  as  much  as  possible  in 
size  but  as  little  as  possible  in  shape  or  polarity.  The  series  n-pentane, 
n-heptane  n-decane.  n-tetradecane  fits  these  specifications  almost  exactly. 

A  plot  of  E  against  the  molar  volume  V  for  these  substances  is  almost 
linear  (n-heptane  falls  slightly  below  the  straight  line  that  passes  through  the 
points  for  n-pentane,  n-decane,  and  n-tetradecane)  with  a  slope  of  0.0348 
kcal/mole/cc.  That  is  to  say.  an  increase  of  1  cc  in  the  molar  volume  that 
leaves  the  overall  shape  and  polarity  unchanged  increases  E  by  0.  0348kcal/ 
mole;  correspondingly,  a  change  of  100  cc  in  V  will  alter  E*p  by  3.48  kcal / 
mole.  This  information  may  also  be  expressed  as  follows: 

(BE^/Bsize)  ,  ,  .  r  0.0348  kcal/mole/cc 

P  shape  polarity 

or,  since  we  shall  use  V  as  a  measure  of  size,  V /L  as  a  measure  of  shape 
and  aaH  as  a  measure  of  polarity, 

(7-1)  (BEp /dV)V/,L  =  0.0348  kcal/mole/cc. 

With  this  information,  we  are  now  in  a  position  to  investigate  the  ef¬ 
fect  of  polarity  upon  Ep. 

8.  The  Dependence  of  Ep  Upon  Polarity. 

The  strategy  this  time  is  to  consider  a  series  of  compounds  that  have 
essentially  the  same  shape  but  which  differ  in  polarity.  They  may  also 
differ  in  size,  in  which  case  we  "correct"  the  observed  Ep's  for  differences 
in  molar  volume  using  equation  (7-1).  Thus  we  could,  for  example,  choose 
among  the  linear  molecules  say  dibutyl  ether  as  a  reference  for  size 
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(V  =  170  cc)  and  say  that  relative  to  this  molecule  the  E_  for  say  ethyl  ace¬ 
tate  (V  =  97.9)  should  be  smaller  by  an  amount  (170  -  98)  (0.0348)  =  2.5 
kcal/mole  if  there  were  no  differences  in  shape  or  polarity  between  the  two. 
But,  although  the  shape  factor  V /L  is  not  too  greatly  different,  there  are, 
however,  significant  differences  in  polarity,  the  MH  index  for  ethyl  acetate 
being  2.4  whereas  that  for  the  ether  is  only  0.8. 

This  difference  in  polarity  tends  to  increase  the  Ep  value  for  ethyl 
acetate  relative  to  that  for  dibutyl  ether.  In  fact,  the  ester's  E  value  is 
greater  than  that  of  the  ether  by  17.7  -  15.9  =  1.8  kcal/mole.  %e  say 
therefore  that  a  change  of  2.4  -  0.  8  -  1.  6  in  the  polarity  factor  increases 
Ep  by  2.5  +  1.8  :  4.  3  kcal/mole,  or  Ep  increases  by  4.  3/1.  6  -  2.7  kcal/ 
mole  per  unit  change  inAAH  (also  measured  in  kcal/mole). 

Actually,  in  reducing  the  data  for  the  straight -chain  hydrocarbons, 
dibutyl  ether,  the  two  acetates,  and  methyl  ethyl  ketone  for  volume  differ¬ 
ences,  it  is  easier  to  calculate  (EpObserved.  0.0348V)  for  each  case.  This 
has  been  done  and  the  values  plotted  against  MH,  the  polarity  factor.  The 
slope  of  the  best  straight  line  through  the  points  which  scatter  a  bit  but  not 
as  much  as  expected  is  2.5.  This  factor  was  then  used  along  with  that  of 
the  previous  section  to  make  a  tentative  estimate  of  the  shape  factor  V /L 
as  illustrated  in  the  following  section.  Using  this  estimate  of  V /L,  one  can 
then  go  back  and  partially  at  least  correct  the  (E^s  -  0.  0348V)  values  used 
above  for  slight  differences  in  the  shape  factor,  obtaining  thereby  a  better 
value  than  2.  5  for  the  polarity  factor  which  can  then  be  used  to  redetermine 
V/L.  etc.  Thus,  by  the  method  of  successive  approximations  we  obtain 
finally  the  value  2.4  for  the  rate  of  change  of  E  with  respect  toAAH,  or 
more  concisely, 

(8-1)  (8E/aAAH)v  -  2.4 


Since  from  Table  6  we  observe  that  variations  of  3-4  kcal/mole  in 
isAH  are  not  unusual,  Ep  may  frequently  vary  from  7.2  -  9.6  kcal/mole  be¬ 
cause  of  the  polarity  factor  alone.  From  equation  (7-1)  one  estimates  that 
it  takes  a  change  in  molar  volume  of  from  7.2/0.0348  =  206  to  9.6/0.0348  = 
276  cc  to  cause  such  a  change  in  Ep.  Such  changes  in  size  are  unusual 
(Table  6),  the  majority  of  the  molar  volumes  fluctuating  by  one-half  to  one- 
third  this  amount,  demonstrating  again  that  the  polarity  factor  is  generally 
two  or  three  times  as  important  as  the  size  factor. 

Now,  as  already  indicated  in  this  section,  we  are  prepared  to  estimate 
the  importance  of  shape,  other  things  being  equal. 

9.  The  Dependence  of  Ep  Upon  Shape. 

The  consideration  of  this  variable  has  been  left  to  last  because,  first 
of  all,  the  size  factor  could  be  obtained  immediately  from  the  straight  chain 
paraffins  without  any  quantitative  information  whatsoever  concerning  either 
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the  shape  or  polarity  factors,  and  secondly  because  prior  information  and 
preliminary  calculations  definitely  indicated  that  of  the  remaining  two  polarity 
was  in  most  instances  the  dominating  factor.  Thus,  whereas  small  vari¬ 
ations  in  the  shape  factor  threw  off  the  initial  &AH  calculation  somewhat, 
even  seemingly  minor  variations  in  AAH  would  make  it  virtually  impossible  to 
draw  reliable  quantitative  conclusions  regarding  the  shape  factor.  That  is  to 
say  the  method  of  successive  approximations  converges  most  rapidly  if  we 
start  with  the  polarity  factor  and  then  proceed  to  the  shape  factor  as  follows. 

First,  however,  it  would  perhaps  be  in  order  at  this  point  to  explain 
why  V /L  was  chosen  as  quantitative  measure  of  the  effect  of  molecular  shape 
upon  Ep. 

Up  to  now  the  molecules  that  have  been  considered- -those  of  sections 
7  and  8- -have  all  been  essentially  rodlike  in  character.  No  mention  has 
been  made  of  the  more  compact  molecules  such  as  ter -butyl  alcohol,  cyclo¬ 
hexane,  or  ortho-xylene  which  in  every  case  (Table  3)  have  a  larger  Ep 
than  the  corresponding  rod -like  molecules. 

This  is  to  be  expected,  for  it  seems  likely  that  the  increase  in  Ep  for 
cyclohexane  over  that  for  n-hexane  arises  from  the  fact  that  in  permitting  a 
molecule  to  diffuse  from  one  equilibrium  lattice  site  to  another  the  polymer 
must  open  up  a  bigger  hole  for  cyclohexane  than  for  n-hexane.  This  is  re¬ 
flected  in  Ep.  That  is,  n-hexane  can  squeeze  through  a  smaller  hole  than 
cyclohexane.  For  this  reason,  an  attempt  has  been  made  in  the  past  to 
correlate  E  values  (actually  the  E^  part  of  Ep)  with  the  dimensions  of  the 
smallest  hole  through  which  the  molecules  can  pass  ("The  Diffusion  of  Some 
Halo-Methanes  in  Polystyrene"  by  G.  S.  Park,  Trans.  Far.  Soc .  46,  684 
(1950)  ).  Although  it  served  to  account  for  some  of  the  experimental  facts, 
this  smallest-hole  theory  was  not  entirely  successful  in  reducing  the  data  to 
a  common  denominator  even  though  relatively  few  compounds  were  studied 
perhaps  in  part  because  polar  effects  which  we  have  found  to  be  so  important 
were  not  considered.  It  seems  likely,  therefore,  that  the  smallest-hole 
theory  really  hasn't  been  given  a  fair  test.  Nevertheless,  we  reject  it  in 
favor  of  the  V /L  index  for  the  following  reasons. 

First  of  all,  it  is  often  difficult  to  clearly  determine  a  "smallest  hole" 
for  a  molecule  using  molecular  models,  for  by  twisting  the  model  around  a 
bit  it  can  be  made  to  slip  through  a  hole  in  a  piece  of  plywood  or  cardboard 
which  very  often  is  significantly  smaller  than  the  smallest  tube  through  which 
the  model  will  pass.  That  is.  the  effect  of  the  diameter  of  your  "smallest 
hole"  depends  upon  just  how  thick  your  hole  is,  and  it  isn't  at  all  obvious 
what  this  thickness  should  be.  Off  hand  one  might  guess  possibly  the  length 
of  two  methyl  groups  (about  4-5  cm  using  the  Fisher -Taylor -Hir schfelder 
models),  but  this  could  only  be  checked  by  constructing  a  whole  set  of  holes 
of  varying  thicknesses  and  diameters  which  could  then  be  used  to  obtain 
several  sets  of  "smallest  hole"  dimensions  to  try  on  the  data.  All  in  ail  a 
fairly  lengthy  project.  Furthermore,  insofar  as  the  following  observations 
are  correct,  it  isn’t  exactly  the  dimension  of  the  "smallest  hole"  we  want 
anyway. 
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For  consider  several  of  the  isomers  of  decane.  According  to  the 
smallest-hole  theory  2,2-dimethyl  octane  and  2, 2,  3,  3-tetramethyl  hexane  will 
have  the  same  energy  of  activation,  for  they  will  pass  through  thd  same 
"smallest  hole"  (and  they  have  essentially  the  same  molar  volumes).  Still 
one  just  sort  of  feels  that  the  tetramethyl  isomer  is  bulkier  than  the  dimethyl 
octane  and  should  experience  greater  difficulty  in  getting  through  polyethylene. 
The  reason  for  this  seems  to  be  that  since  polyethylene  is  a  linear  polymer, 
one  naturally  expects  that  random  thermal  motion  of  the  polymer  segments 
will  open  up  rod  shaped  holes  more  often  than  spherically  shaped  ones  of 
equal  volume.  There  are  perhaps  more  elegant  ways  of  putting  this.  Thus, 
in  comparing  n-decane,  2,2-dimethyl  octane,  and  2, 2,  3,  3-tetramethyl  hexane, 
we  could  say  that  the  standard  free  energy  of  formation  in  polyethylene  of 
banana- shaped  holes  is  less  than  that  for  pear-shaped  holes  which  in  turn  is 
less  than  that  for  small  watermel  on- shaped  holes,  the  volume  of  the  holes 
being  the  same  in  each  case,  namely  194/6x10^3  cm^  (Table  6). 

In  any  event,  whether  we  use  every-day  language  or  that  of  the  Abso¬ 
lute  Reaction  Rate  Theory,  all  we  are  suggesting  at  this  point  is  that  poly¬ 
ethylene  will  in  fact  distinguish  between  the  two  isomers  of  n-decane  mention¬ 
ed  above.  Since  the  smallest-hole  theory  does  not  predict  this,  we  reject  it 
for  the  V/L  parameter  which  does.  Furthermore,  L  is  easily  determined  by 
simple  measurement  although  here  too  there  is  ambiguity  in  some  cases. 
Generally  what  has  been  done  is  to  stretch  the  model  out  straight  and  take 
the  longest  linear  dimension  for  L.  (We  shall  encounter  an  interesting  re¬ 
finement  upon  this,  however,  whan  we  consider  para-derivatives  of  benzene.) 

We  might  summarize  the  argument  leading  up  to  the  choice  of  V  /L  by 
saying  that  this  parameter  gives  us  a  better  measure  of  over-all  bulkiness, 
the  significant  factor,  than  does  the  diameter  of  the  smallest  hole  through 
which  the  molecule  can  pass. 

For  the  straight  chain  hydrocarbons  V/L  is  constant  at  about  12.  The 
presence  of  an  oxygen  atom  in  an  essentially  linear  molecule  decreases  the 
parameter  by  about  a  unit  or  so.  Table  6  shows  that  the  value  of  V/L  most 
often  falls  between  9  and  12  with  water  and  formic  acid  (5.2  and  6.7  re¬ 
spectively)  and  ter-butyl  alcohol,  cyclohexane,  and  ortho-xylene  (15.9,  15.6, 
and  13.  9  respectively)  at  the  extremes. 

To  show  the  essence  of  the  method  used  to  determine  the  dependence 
of  Ep  upon  V/L,  we  compare  in  Table  7  the  parameters  for  ethyl  acetate  and 
cyclohexane  as  given  in  Tables  3  and  6. 

Table  7 


A  Comparison  of  Ethyl  Acetate  and  Cyclohexane. 


Substance 

Vcc 

AAH 

V/L 

0.0348V 

2.4A&H 

E„ 

cc  /cm 

cc 

kcal/moi  kcal/mol 

Ethyl  Acetate 

97.9 

2.4 

10.6 

3.4 

5.  8 

17.  7 

Cyclohexane 

81.4 

0.6 

15.  6 

2.  8 

1.4 

16.2 
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Since  the  ethyl  acetate  molecule  is  both  larger  (column  2)  and  more 
polar  (column  3)  than  the  cyclohexane  molecule,  we  expect  that  (equations 
7-1  and  8-1)  owing  to  these  factors  its  activation  energy  will  be  greater  by 
(3. 4+5.  8) -(2.  8-<-l .  4)  =  5.  0  kcal/mole  (columns  5  and  6)  whereas  in  fact  (last 
column)  the  difference  is  only  17.7-16.2-1.5  kcal/mole  in  this  direction. 

We  say  then  that  the  E_  for  cyclohexane  has  been  raised  relative  to  that  for 
the  acetate  ester  by  5. 0-1.  5-3.  5  kcal/mole  because  it  is  bulkier  (column  4) 
by  15.6-10.6  =  5.0  units.  That  is,  in  this  case  the  rate  of  change  of  Ep  with 
respect  to  V /L  is  3.5/5  =  0.70.  F 

A  number  of  substances  (ethyl  acetate,  methyl  ethyl  ketone,  dibutyl 
ether,  acetone,  the  paraffins,  and  cyclohexane)  have  been  compared  in  es¬ 
sentially  this  manner  by  plotting  (Ep-0.0348V-2.4AAH)  against  V /L.  The  best 
straight  line  through  the  points  has  a  slope  of  0.  75,  so  that  we  have  finally 
to  go  along  with  equations  (7-1)  and  (8-1)  the  following  value  for  the  shape 
factor : 

(9-1)  OEp/dV/L)V(AAH  =  0.75 

As  already  observed,  the  parameter  V  /L  generally  falls  between  9 
and  12.  Hence  the  contribution  the  shape  factor  makes  to  Ep  may  often  by 
expected  to  vary  by  about  3(0.  75)  =2.  2  kcal/mole  though  seldom  by  twice  this 
amount.  Comparing  this  with  the  results  of  the  previous  sections,  we  con¬ 
clude,  therefore,  that  generally  speaking  the  polarity  factor  is  most  important 
followed  by  size  and  then  shape.  These  three  factors  are  synthesized  into  a 
single  expression  in  the  following  section. 


10.  The  Quantitative  Dependence  of  Ep  Upon  Size,  Shape,  and  Polarity. 

To  summarize  the  argument  up  to  this  point  concerning  the  structure  of 
Ep,  it  was  suggested  in  section  6  that  the  observed  activation  energy  for 
permeation  should  depend  upon  three  factors  relating  to  the  properties  of  the 
permeant  molecule,  namely  its  size,  shape,  and  polarity.  hi  principle  then 
variations  in  the  observed  Ep  values  are  decomposable  into  three  parts  as 
indicated  below. 

(10-la)  dEp  =  OEp/Ssize)shapef  d(size)  +  (0Ep/3shape)size^  d(shape) 

polarity  polarity 

-+•  (8Ep/SP°larity)size(  d(  polarity). 

shape 

If  we  take  the  molar  volume  V  (the  gram  molecular  weight  divided  by 
the  density)  as  a  measure  of  size,  the  parameter  V/L  (L,  =  the  length  of  the 
stretched-out  Fisher- Taylor -Hir schfelder  molecular  models)  as  an  indication 
of  over-all  shape,  and  AAH  (defined  in  equation  5-3)  as  an  index  of  polarity, 
equation  (10- la)  becomes 


WADC  TR  53-133  Pt  2 


24 


(10- lb) 


dEp  =  OEp/aV)v/L,A&HdV  +  (9Ep/3V/L)v>MHd(V/L) 

+  (8Ep^AAH)y>  v/Ld(^AH), 

In  the  preceding  sections  (7,  8,  and  9),  it  was  shown  how  the  data 
of  this  project  could  be  manipulated  to  yield  values  for  the  differential  coef¬ 
ficients  appearing  in  equation  (10- lb).  Using  these  results  (equations  7-1, 
8-1,  and  9-1),  we  obtain  finally 

(10-  lc)  dEp  =  0.  0348dV  •+  0.75d(V/L)  4-  2.4d£AH). 

The  problem  now  is  to  integrate  this  equation.  Proceeding  directly  as 
indicated,  we  get 

(10-2a)  Ep  =  0.0348V -i-  0.  75(V/L)  +  2.4(AAH)  -t-  constant 

where  the  constant  of  integration  remains  to  be  evaluated.  Note,  however, 
that  this  can  not  be  done  simply  by  arguing  that  Ep->0  as  V  -»  0  since  V/L 
and  A  AH  do  not  necessarily  approach  zero  with  V  even  if  Ep  should,  and,  in 
fact,  the  evidence  is  against  even  this  too. 

One  could,  of  course,  proceed  empirically  at  this  point,  substituting 
into  (10-2a)  from  Tables  3  and  6.  If  concordant  results  were  obtained,  this 
procedure  would  not  only  yield  a  reliable  value  for  the  integration  constant, 
but  also,  it  would  seem,  should  serve  as  a  check  on  the  assumption  inherent 
in  the  derivation  of  (10-2a)  that  the  differential  coefficients  in  (10- lb)  are 
independent  of  the  values  for  V,  V  /L,  and  A  AH.  Actually,  this  seems  un¬ 
likely.  Thus,  for  example,  one  expects  that  the  rate  of  increase  in  the 
energy  required  to  form  a  hole  with  respect  to  increasing  V  or  V/L  should 
itself  increase  with  increasing  AAH.  That  is  to  say,  one  suspects  that  actu¬ 
ally  the  differential  coefficients  (3Ep/3V)y /Lt aaH  an(^  OEp/SV 

are  increasing  functions  of  AAH  rather  than  constants  with  the  values  0.0348 
and  0.75  respectively.  This  does  not  look  too  good  for  the  theory.  Fortu¬ 
nately  there  is  a  rather  elegant  way  out  of  all  this. 

First  we  observe  that  although  E  increases  with  say  V,  Ep  is  not 
proportional  to  V  as  can  be  seen  from^the  fact  that  Ep  increases  from  12.8 
to  only  15.6  in  going  from  pentane  to  decane;  that  is,  E_(2V,  V /L,AAH)  ^ 
2Ep(V,  V/L, AAH).  If,  however,  at  the  same  time  V  is  doubled  V/L  and  AAH 
are  also  doubled,  and  if  then  Ep  doubles,  or,  more  generally,  if 

( 10-2b)  Ep(AV,XV/L,XAAH)  =  XEp(V,  V/L,AA  H) 

where  X  is  any  arbitrary  positive  number,  it  follows  immediately  from 
Euler's  theorem  for  homogeneous  functions  that 

<  1 0 -2c)  Ep  =  (8Ep/av)v/L>AfikHv  +  (aEp/3V/L)yiAAH(V/L) 

+  (3Ep/aAAH)yy/L(AAH). 
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The  discussion  leading  to  equation  (10-2c)  does  not  require  that  the 
differential  coefficients  be  constant  over  the  whole  range  of  possible  V,  V/L, 
andAAH  values.  Nevertheless,  within  the  range  of  experimental  interest 
the  coefficients  do  indeed  appear  to  be  nearly  constant  in  value  as  indicated 
by  equations  (7-1),  (8-1)  and  (9-1).  Substituting  from  these  into  (10-2c),  we 
obtain  finally 

(10~2e)  Ep  =  0.0348V  +  0.75V/L  +  2.4&AH 

where  the  integration  constant  of  equation  (10-2a)  appears  with  a  "theoreti¬ 
cally"  assigned  value  of  0.  In  the  following  section  we  begin  a  comparison 
of  equation  (10~2e)  with  experiment. 


11.  The  Calculation  of  E^  Values  for  the  Aliphatic  Hydrocarbons, 

Ethers,  Esters,  and  Ketones. 

Using  equation  (10-2e)  and  the  data  of  Table  6,  Ep  values  for  straight- 
chain  compounds  with  polarities  ranging  from  0.0  to  3.6  are  calculated  and 
compared  with  experiment  as  shown  in  Table  8. 

Table  8 

Theory  and  Experiment  Compared  for  Straight -Chain  Compounds 
With  Polarities  Between  0  and  3. 


Substance 

0. 0348V 

0.75V  /L 

2.4aaH  E 

calc. 

E  6XP' 

Diff. 

cc 

cc  /cm 

kcal/mole 

kcal/mol^  kcal/mol 

n-Pentane 

4.0 

9.0 

0.0 

13.  0 

12.  8 

+0.2 

n- Heptane 

5.  0 

9-0 

0.  0 

14.  0 

13.2 

-+0.  8 

n- Decane 

6.  7 

9.  0 

0.  0 

15.  7 

15.  6 

+•0.  1 

n-Tetradecane 

9.0 

9.0 

0.  0 

18.  0 

18.  0 

+0.  0 

Dibutyl  ether 

5:9 

8.  5 

1.9 

16.  3 

16.  0 

+0.  3 

Cyclohexane 

3.7 

11.7 

1.4 

16.  8 

16.2 

*0.  6 

Ethyl  Acetate 

3.  4 

8.  0 

6.  7 

18.  1 

17.  7 

+0.4 

ME  Ketone 

3.  1 

8.  5 

6.  7 

18.  3 

18.0 

4.0.  3 

Acetone 

2.  6 

8.  5 

7.  2 

18.  3 

18.  7 

-0.4 

Amyl  Acetate 

5.  1 

8.4 

6.  0 

19.  5 

19.4 

4-0.  1 

There  are  several  things  of  interest  in  Table  8.  First  of  all, 

a.  The  hypothesis  of  a  value  zero  for  the  integration  constant 
of  equation  (10 -2a)  is  confirmed.  In  fact,  since  as  recorded  in  the  last 
column  the  calculated  value  is  usually  slightly  greater  than  the  observed  value, 

b.  The  overall  fit  of  the  data  by  equation  (10-2e)  could  undoubted¬ 
ly  be  improved  a  bit  by  a  slight  readjustment  in  the  three  parameters  0.0348, 
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0.75,  and  2.5.  Nevertheless, 


c.  The  three  parameter  equation  generally  reproduces  the  ten 
observed  values  to  within  2-3%.  Hence  it  is  felt  that 


d.  Equation  (10-2e)  can  certainly  reliably  predict  the  order  of 
magnitude  of  activation  energies  for  the  aliphatic  hydrocarbons,  ether s,  esters, 
and  ketones.  By  the  order  of  magnitude  of  Ep  we  mean  that  the  correspond¬ 
ing  P  values  (equation  1-1)  will  not  be  in  error  by  more  than  a  factor  of  ten. 
From  equation  (1-2)  it  follows  that  at  room  temperature  this  amount  of 
leeway  requires  that  E  be  correct  to  within  1.4  kcal/mole.  For  an  accuracy 
within  a  factor  of  two,  ^the  Ep  calculation  must  not  be  in  error  by  more  than 
0.4  kcal/mole.  At  higher  temperatures  less  accuracy  is  required.  Also,  the 
method  for  estimating  Pq  values  (to  be  presented  later)  automatically  compen¬ 
sates  in  part  for  errors  in  Ep.  Thus  agreement  between  theory  and  experiment 
better  than  that  suggested  above  would  not  be  surprising.  Finally, 


e.  The  structure  of  E  is  nicely  revealed  by  the  first  three 
columns  of  figures.  For  example,^  the  Master  Plot  shows  that  tetradecane 
and  ME  ketone  exhibit  similar  permeability  behavior  but,  as  Table  8  shows, 
for  rather  different  reasons:  the  size  and  polarity  factors,  operating  in  op¬ 
posite  directions,  tend  to  cancel. 


12. 


The  Calculation  of  E^  Values  for  the  Aromatic  Hydrocarbons,  the 
Substituted  Aromatic  Hydrocarbons,  and  Dipentene. 


In  this  section  we  consider  to  what  extent  our  formula  for  E^  is  ap¬ 
plicable  to  the  aromatic  compounds  and  the  cyclic  compound  dipentene. 


c  sic 

The  values  listed  for  E  as  E  in  column  3  of  Table  9  are  calcu- 

p  p 

lated  using  equation  (10-2e)  and  the  data  of  Table  6.  These  are  followed  by 
the  experimental  values  in  column  4  and  then  a  column  headed  EpCorr  to  be 
explained  later. 


Table  9 


Theoretical  and  Experimental  Ep  Values  for  Some 

Aromatic  Compounds 


Compound  Ref.  Cpd. 

for  A&H  Calc. 

E  calc 
P 

E  exp 
P 

£  corr 
P 

E  corI  E  eXp 

P  P 

Benzene 

Hexane 

13.  1 

13.  3 

13.  1 

-0.2 

o- Xylene 

Octane 

16.  0 

15.  6 

16.0 

r0. 4 

p- Xylene 

Octane 

14.7 

13.4 

13.7 

«-0.  3 

Chlorobenzene 

Heptane 

15.6 

14.2 

14.  3 

+  0.  1 

p-Chlorotoluene 

Octane 

17.2 

15.  8 

15.4 

-0.4 
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The  calculated  value  for  benzene  is  probably  within  experimental  error, 
that  for  o-xylene  a  bit  high,  and  p-xylene  is  considerably  off  (1,3  kcal/mole). 
Chlorobenzene  and  p-chlorotoluene  are  also  each  too  high  by  1.4  kcal/mole. 

These  discrepancies  may  merely  reflect  second  order  structural  ef¬ 
fects  that  are  inadequately  accounted  for  by  a  single  shape  factor  such  as 
V /L.  The  structures  of  these  molecules  are  shown  below  in  Figure  4. 


Benzene 


o-Xylene  p-Xylene  Chlorobenzene 


Fig.  4  Structures  of  Aromatic  Compounds 


p-Chlorotoluene 


Now  presumably  once  a  hole  has  been  opened  up  large  enough  to  ac¬ 
comodate  a  diffusing  benzene  molecule,  it  makes  little  difference  if  we  then 
attach  methyl  groups  fore  and  aft  as  in  p-xylene  since  this  is  probably  "dead 
air  space"  anyway;  p-xylene  is  so  to  speak  "streamlined". 

Similarly,  the  chlorine  atom  which  tags  along  with  the  benzene  nucleus 
in  chlorobenzene  probably  is  acting  in  a  partial  vacuum  so  far  as  its  polarity 
effects  are  concerned. 

If  we  apply  a  correction  of  0.5  kcal/mole  for  each  methyl  group  in  p- 
xylene  and  1.3  kcal/mole  for  the  chlorine  atom  in  chlorobenzene,  we  obtain 
so-called  corrected  values  as  listed  in  the  fifth  column  of  Table  9.  Applying 
both  of  these  corrections  to  the  last  molecule,  p-chlorotoluene  also  brings 
its  calculated  value  into  reasonable  agreement  with  experiment. 

With  respect  to  these  calculations,  we  observe  that 

a.  Although  they  are  undoubtably  based  upon  an  over-simplified 

model,  and 


b.  Although  they  are  highly  empirical, 

c.  It  is  nevertheless  a  fact  that  whereas  the  first  two  values  in 
Table  9  agree  about  as  usual  with  experiment,  the  last  three  are  off  by  more 
than  two  or  three  times  the  usual  amount. 
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d.  Furthermore,  it  might  be  mentioned  that  the  correction  sug¬ 
gested  by  the  p-xylene  -  benzene  comparison  did  not  at  first  seem  to  account 
for  results  obtained  on  chlorobenzene  and  p-chlorotoluene.  A  check  of  all 
the  calculations  involved  finally  revealed  an  error  in  the  1  mm  column  of 
Chart  1  on  page  16  of  Dreisbach. 

e.  Finally,  this  "streamlined"  effect  also  seems  to  be  operative 
in  the  case  of  the  dipentene  molecule  pictured  below  in  Figure  5. 


/ 

C 

\\ 


Fig.  5.  Structure  of  Dipentene 

Treating  this  compound  as  a  hydrocarbon  (A AH  =  0),  we  calculate  that  E  = 
0.0348(162)  +  0.75(15.  1)  =  16.9  kcal/mole  as  compared  to  the  observed  value 
of  16.0  kcal/mole. 

In  general,  then,  we  expect  that  para- substituted  compounds  will  have 
Ep  values  about  1  kcal/mole  less  than  that  calculated  using  equation  (10-2e). 

The  following  section  .gives  some  calculations  on  the  more  highly  polar, 
strongly  hydrogen-bonded  alcohols  and  acids  whose  &AH  values  axe  greater 
than  3.  0. 

13.  Calculated  and  Experimental  E^  Values  for  the  Alcohols  and  Acids. 

Table  10 


Substance 

0.  0348Y 

cc 

0.  75V /L 
cc/cm 

2.4MH 
k/c  al  /mol 

E  calc  E  exp  Difference 
k /cal /mol  1  /cal  /mol 

n-Propanol 

2.  6 

7.  2 

14.  9 

24.  7 

22.  1 

+  2.6 

n- Butanol 

3.2 

7.  5 

14.  0 

24.  7 

22.4 

+  2.  3 

sec  -  Butanol 

3.2 

8.9 

12.  9 

25.0 

23.0 

2.  0 

ter -Butanol 

3.  3 

11.9 

11.  5 

26.  7 

26.2 

+  0.  5 

n-Octanol 

5.  5 

8.  1 

12.2 

25.  8 

22.5 

+*  3.  3 

Phenol 

3.0 

7.  7 

12.  5 

23.2 

21.2 

'2.0 

Acetic  Acid 

2.0 

6.  8 

14.  6 

23.4 

17.  8 

-  5.  6 

Butyric  Acid 

3.2 

6.  5 

17.  7 

27.4 

18.  3 

+"9.  1 
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The  explanation  for  the  relatively  large  discrepancies  between  the 
calculated  (equation  10-2e)  and  experimental  values  is  probably  this. 

Suppose  these  substances  dissolve  in  polyethylene  as  dimers  rather 
than  single  monomeric  molecules.  (Additional  evidence  for  believing  this 
will  be  presented  later.)  Then  when  it  comes  to  diffusing  from  one  equi¬ 
librium  position  to  another  there  exist  at  least  two  possibilities: 

(1)  The  dimer  splits  apart  as  the  transition  state  for 
diffusion  is  approached,  or 

(2)  The  dimer  as  a  whole  diffuses  from  one  equilibrium 
position  to  another. 

Now  presumably  all  polar  substances  are  to  some  extent  associated 
(dimerized)  in  polyethylene  which  so  to  speak  is  a  poorly  ionizing- -or .better, 
poorly  dissociating- - solvent;  i.  e.  ,  it  has  no  great  tendency  to  solvate 
(dissolve)  the  solute.  Thus  solute- solute  contacts  as  pictured  in  Fig.  6  are 
energetically  favored  over  solute- solvent  (polyethylene)  contacts,  and  at  abso¬ 
lute  zero  (no  random  thermal  motion)  dimerization  would  be  complete.  (In 
fact,  at  this  temperature  even  the  dimers,  which  undoubtably  still  exhibit  a 
slight  residual  polarity,  would  tend  to  dimerize,  etc.  until  eventually  a  pure 
phase  begins  to  separate  out,  merely  another  way  of  saying  that  at  T:0  the 
solubility  is  zero.  ) 


Figure  6 

Mutual  Attraction  of  Dipoles  in  Nitrobenzene 
(G.  Briegleb,  Zwischenmolekulare  Kraftz  and  Molekulstruktur  Pg.  55) 


But  at  any  finite  temperature  random  thermal  motion  (the  entropy  ef¬ 
fect)  will  tend  to  break  up  the  dimers  until  the  free  energy  (E-TS)  is  a  mini¬ 
mum.  For  the  slightly  polar  ethers  this  thermal  breakdown  is  probably 
nearly  100%  complete  at  room  temperature.  There  is  no  reason  in  this  case 
for  suspecting  that  the  diffusing  species  is  anything  other  than  the  simple 
monomeric  molecule. 

With  a  slightly  more  polar  substance  such  as  acetone,  however,  there 
may  be  some  association,  both  in  the  pure  liquid  and  in  the  liquid  dissolved 
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in  polyethylene.  On  the  other  hand,  the  gas  in  this  case  is  not  associated  so 
that  the  full  polarity  of  acetone  shows  up  in  <\/\H  as  it  will  too  in  E  if  it  is 
the  monomer  that  diffuses  rather  than  the  dimer;  that  is,  if  it  is  thermo¬ 


dynamically  easier  to  break  the  not-too-strong  intermolecular  bonds  between 
two  acetone  molecules  and  push  through  a  single  molecule  of  acetone  rather 
than  conserve  at  the  first  point  by  leaving  such  dimers  as  there  are  alone 
and  then  trying  to  push  through  this  bulkier  unit  later  on.  In  cases  such  as 
this  the  reasoning  leading  up  to  the  introduction  of  the  polarity  factor  &AH  is 
valid  and  we  expect  equation  (10-2e)  to  hold. 


Still,  there  must  eventually  come  a  point  where  it  is  thermodynamically 
uneconomical  to  sever  the  dimer  before  diffusion.  Tetradecane,  for  example, 
undoubtably  diffuses  as  just  that:  a  14  carbon  hydrocarbon,  rather  than  as 
two  "heptane"  fragments.  This,  of  course,  is  an  extreme  case.  Indeed,  the 
"intermolecular  force"  between  two  seven  carbon  tetradecane  fragments  is 
essentially  a  strong,  localized,  honest-to-goodness  100%  chemical  bond.  This 
bond  (along  with  all  the  other  bonds  in  the  molecule)  is  so  strong  that  even 
the  possibility  of  great  randomness  in  the  gaseous  state  fails  to  cause  the 
molecule  to  dissociate  upon  vaporization.  Whatever  "polarity"  we  might  wish 
to  ascribe  to  a  tetradecane  fragment  is  thus  completely  hidden  from  us  at  all 
times.  At  no  point  is  there  ever  any  excuse  for  including  the  great  "polarity" 
of  any  tetradecane  fragment  in  the  AAH  for  this  substance. 


With  the  alcohols  and  the  acids,  however,  this  is  not  the  case.  The 
intermolecular  bonds  here  are,  relative  to  the  hydrocarbons,  ethers,  esters, 
and  ketones  very  strong,  so  strong,  in  fact,  that  they  earn  a  special  name: 
hydrogen  bonds 

\c-C  C  -fc, 

- ->  C  ' 


Figure  7 

Hydrogen  Bond  in  Acids 


Figure  8 

Hydrogen  Bond  in  Alcohols 
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But  the  hydrogen  bonds  in  the  alcohols  and  acids  are  not  strong  e- 
nough  to  maintain  a  degree  of  dimerization  in  the  gaseous  state  anything  like 
that  in  the  liquid  state.  Hence  pretty  much  the  full  polarity  of  these 
substances  shows  up  in  the  heat  of  vaporization  and,  therefore,  inAAH.  When 
it  comes  to  diffusing  in  the  polymer,  however,  thisAAH  factor  is  weighted 
rather  heavily  (the  factor  is  2.4),  and  in  these  alcohols  and  acids  it  is  ap¬ 
parently  thermodynamically  easier  to  diffuse  the  dimer  rather  than  the  more 
polar  although  smaller  monomer. 

If  we  knew  the  shape  of  the  dimer  (Figure  6)  and  its  residual  polarity, 
it  should  be  possible  to  calculate  by  equation  (10-2e)  a  value  for  Ep  which 
then  could  be  compared  with  that  for  the  simple  molecule  itself.  Other  things 
being  equal  (namely  Pq),  the  observed  Ep  would  be  the  lower  of  the  two.  Un¬ 
fortunately  neither  the  shape  nor  the  polarity  of  the  dimer  are  known. 

With  these  remarks  in  mind,  we  turn  back  to  a  consideration  of  Table 

10. 


Among  the  alcohols  the  calculated  value  is  on  the  average  about  2.  1 
kcal/mole  too  large  except  in  the  case  of  ter -butanol  which,  owing  to  both  its 
great  bulkiness  and  reduced  polarity,  probably  still  diffuses  as  the  dimer. 

The  aromatic  alcohol  phenol  behaves  like  the  aliphatic  alcohols  in  so  far  as 
Ep  is  concerned. 

The  two  acids  in  Table  10  deviate  from  "ideal"  behavior  even  more  than 
the  alcohols,  reflecting,  no  doubt,  the  well  established  fact  that  hydroge’n 
bonding  in  the  carboxylic  acids  is  stronger  than  in  the  corresponding  alcohols 
(Table  4).  Indeed,  the  hydrogen  bonding  in  acetic  acid  is  known  to  carry 
over  to  some  extent  to  the  gaseous  state,  hiding,  therefore,  some  of  its  polar¬ 
ity  as  indicated  byA/\H.  It  might  well  be  that  a  study  of  gas  imperfections 
and  the  temperature  variation  of  the  second  virial  coefficient  would  reveal  the 
reasons  for  the  discrepancies  between  acetic  and  butyric  acid  recorded  in  the 
last  column  of  Table  10. 

Table  11 


Comparison  of  Some  Calculated  and  Observed  Ep  Values 


Material 

Calc.  Ep  k/cal/mol 

Obs.  Ep  k/cal/mol 

Diacetone  Alcohol 

26 

21.9 

Allyl  Alcohol 

23 

20.2 

Methyl  Alcohol 

26.  7 

17.7 

W  ater 

23.  5 

18.4 

Formic  Acid 

20.  2 

16.7 

Benzaldehyde 

17.9 

19.  3 

B  utyr  aldehyde 

16.  7 

19.0 
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In  diacetone  alcohol  the  difference  (26.0-21.9=4.1)  is  larger  than  in 
the  other  alcohols  because  the  ketone  oxygen  may  be  partially  screened  in  the 
dimer;  allyl  alcohol  is  not  far  off  from  the  other  alcohols,  and  the  great 
decrease  (26.7-17.7=9.0)  found  for  methyl  alcohol  is  probably  explained  by 
the  tendency  of  this  relatively  small  molecule  to  trimerize,  etc.*  Water  also 
probably  tends  to  form  higher  polymers  although  there  are  no  methyl  groups 
to  screen  the  oxygens  in  this  case,  and  the  decrease  in  Ep  is  only  5.  1  kcal/ 
mole.  Formic  acid  continues  the  trend  already  observed  with  the  other  two 
acids.  Finally,  the  calculated  values  for  the  relatively  non-polar  aldehydes 
are  strangely  too  low  by  1.  3-1.4  kcal/mole. 


14. 


The  Calculation  of  E 

- p 


A  Summary. 


From  the  above  and  earlier  remarks,  it  is  abundantly  evident  that  the 
penetration  of  polyethylene  by  organic  solvents  is  a  very  selective  process, 
one  which  not  infrequently  is  significantly  dependent  upon  seemingly  unim¬ 
portant  structural  features.  Nevertheless,  it  now  seems  possible  to  predict 
quite  good  activation  energies  for  the  aliphatic  hydrocarbons,  ethers,  esters, 
and  ketones  with  only  a  three  parameter  equation,  and  with  the  equation 


(14-1)  Ep  =  0.0348V  +  0.75V/L  +  2.44AH  -  C 

reasonably  reliable  predictions  can  be  made  for  the  substituted  aromatics, 
aldehydes,  and  the  higher  alcohols  provided  the  following  values  are  used  for 
C. 


Table  12 

Values  for  the  Parameter  in  Equation  (14-1) 


Class  of  Compounds  Value  to  Use  for  C 


Straight  Chain  Hydrocarbons,  sat'd  0 

Branched  Hydrocarbons,  sat’d  0 

Unsaturated  Hydrocarbons  0 

Cyclic  Hydrocarbons  0 

Ethers  0 

Esters  0 

Ketones  0 

Aromatic  Hydrocarbons  0 

par  a- substituted  Compounds  1 

Alcohols  2.  3 

Aldehyde  s  -1.3 


*  Propanol  may  be  doing  this  a  bit  too. 


WADC  TR  53-133  Pt  2 


33 


If  information  on  one  or  two  more  acids  were  available  equation  (14-1) 
could  probably  be  extended  to  this  class  of  compounds,  which  has  not  been 
included  in  the  table  below. 


Table  13 

Comparison  of  Calculated  and  Observed  E 


Substance 

£  calc 

kc  al/mol 

Eexp 

kc  al/mol 

£  calc  £  exp 

P  '  P 

n-Pentane 

13.  0 

12. 8 

+  0.2 

n- Heptane 

14.  0 

13.2 

+  0.8 

n-Decane 

15.  7 

15.  6 

+  0.  1 

n-Tetradecane 

18.  0 

18.  0 

0.  0 

Cyclohexane 

16.8 

16.2 

t  0.6 

Dibutyl  Ether 

16.  3 

16.  0 

+  0.  3 

Ethyl  Acetate 

18.  1 

17.  7 

r  0.4 

Amyl  Acetate 

19.  5 

19.4 

+  0.  1 

Acetone 

18.  3 

18.7 

-0.4 

Methyl  Ethyl  Ketone 

18.  3 

18.0 

+  0.  3 

Benzene 

13.  3 

13.  1 

+  0.  2 

o- Xylene 

16.  0 

15.  6 

+■0.4 

p- Xylene 

13.  6 

13.4 

+0.  2 

Chlorobenzene 

14.  3 

14.  1 

+0.2 

p  -  Chlor  otoluene 

16.  2 

15.  8 

+  0.  4 

Dipentene 

15.  9 

16.  0 

-0.  1 

Butyr  aldehyde 

19.  0 

19.0 

o.  o 

Benz  aldehyde 

19.2 

19.  3 

-0.  1 

n -Propanol 

22.4 

22.  1 

*  0.  3 

n-Butanol 

22.4 

22.4 

0.  0 

sec  -Butanol 

22.7 

23.  0 

-0.  3 

ter -Butanol 

26.  7 

26.2 

+  0.  5 

n-Octanol 

23.  5 

22.  5 

+  1.0 

Phenol 

20.9 

21.2 

-0.  3 
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The  average  error  in  the  calculated  value  for  the  activation  energies 
of  the  two  dozen  compounds  listed  in  Table  13  is  0.2^  kcal/mole  although 
this  value  should  not  be  taken  too  seriously  inasmuch  as  the  table  was  con¬ 
structed  with  these  compounds  very  much  in  mind.  In  practice  it  might  not 
be  too  optimistic  to  expect  that  application  of  equation  (14-1)  to  similar 
substances  will  yield  E  values  reliable  to  within  twice  this  average  error  or 
0.  5  kcal/mole  on  the  average.  As  already  mentioned,  to  know  P  itself  with¬ 
in  a  factor  of  ten,  it  would  be  well  to  be  off  by  not  more  than  1.4  kcal/mole 
on  Ep. 

15.  Interlude. 


The  discussion  up  to  now  has  been  concerned  almost  entirely  with  the 
parameter  E_  of  equation  (l-3);nothing  yet  has  been  said  about  PQ  although 
from  the  last  column  of  table  3  it  appears  that  this  parameter  varies  much 
more  from  compound  to  compound  than  does  Ep.  This  observation,  however, 
overlooks  the  exponential  dependence  of  P  upon  E  .  As  shown  in  equation  (1-3) 
it  is  the  logarithm  of  PQ  rather  than  PQ  itself  thlt  should  be  compared  with 
E  .  Indeed,  we  have  observed  already  that  at  room  temperature  a  change 
oF  1.4  kcal/mole  in  E  alters  P  tenfold.  A  spread  in  E  from  26.2  (ter- 
butyl  alcohol)  to  12.8  ^n-pentane)  represents,  therefore,  Ipproximately  a 
10*®  fold  variation  in  the  exponential  factor  of  equation  (1-1)  whereas  the 
factor  PD  varies  only  somewhat  less  than  10^  fold.  (The  variations  of  the  two 
factors  work  in  opposition  so  that  P  itself  at  room  temperature  shows  a 
spread  of  about  0.26  (ter -butyl  alcohol)  to  526  (n-pentane);  about  a  2000  fold 
change . ) 

This  then  was  one  of  the  principle  reasons  for  discussing  Ep  first.  It 
is  1000  times  more  important  than  PQ.  With  information  on  Ep,  furthermore, 
a  single  measurement  of  P  at  some  convenient  temperature  will  enable  one 
to  predict  P-Factors  from  say  -40°F  to  212°F,  for  knowing  P  at  some  T  and 
also  Ep,  one  can  solve  either  equation  (l-2a)  or  (1-3)  for  PQ.  Of  course  PQ 
if  known  along  with  one  P -Factor  would  serve  to  determine  Ep  rather  than 
vice-versa.  But  the  factors  that  determine  PD  are  not  as  familiar  nor  as 
well  understood  as  the  effect  of  size,  shape  and  polarity  upon  E  . 


In  short,  the  problem  of  E  is  one  of  energetics.  It  resolves  itself 
into  concrete  geometrical  factors1^  and  familiar  inter  molecular  interactions. 

The  problem  of  PQ,  on  the  other  hand,  is  essentially  one  of  entropy.  It 
generally  involves  one  in  numerous  difficulties  about  which  unfortunately  there 
is  little  concrete  or  familiar.  To  start  the  thermodynamical  discussion  of 
Po>  we  b  egin  by  considering  the  absorption  tests,  which  are  also  of  interest 
from  the  point  of  view  of  the  activation  energy  E  . 


16.  The  Master  Solubility  Data. 


These  data  have  already  been  reported  by  Plax  as  the  grams  of  pene¬ 
trant  absorbed  per  gram  of  polyethylene,  herein  called  Q.  For  the  slightly 
soluble  solvents  Q  is  almost  identical  with  the  weight  fraction  of  the  pene- 
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trant  in  the  (slightly)  swollen  polymer.  What  we  require  in  all  cases,  how¬ 
ever,  is  the  volume  fraction  (<^>).  If  we  let  p  stand  for  the  polymer  and  2 
for  the  penetrant,  the  two  volume  fractions  are,  respectively, 


(16-1) 


^p  =  V?/(npVp  +  n2V2)  and 


<j)2  =  («2vf) /<npvp  +  n2v2> 


23 

where  the  n's  are  the  mole  numbers  (n.  6.  10  =number  of  particles)  and  the 
v°'s  represent  the  molar  volumes  of  the  pure  components.  Naturally <$>  p  + 
<1>2  =  1- 

The  conversion  from  Q  to  say^p  is  given  quite  simply  by 
(16-2)  «{)  =  1/(1  +  DpQ/D2). 


For  the  density  D  of  polyethylene  we  have  used  the  following  values: 

0.935,  0.922,  0.<#0,  and  0.898  at  32°F,  70°F,  100°F,  and  130°F  respective 
ly.  The  densities  of  the  penetrants  were  calculated  from  data  given  in  the 
International  Critical  Tables. 

It  might  be  observed  at  this  point  that  equations  (16-1)  and  (16-2)  as¬ 
sume  the  volume  of  the  swollen  polymer  can  be  calculated  directly  from  the 
volumes  of  the  pure  starting  materials.  Experimental  evidence  bearing  on 
this  will  be  presented  later. 

In  the  table  below  the  substances  are  listed  approximately  in  order  of 
increasing  solubility  in  polyethylene  at  70°F. 
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,  Table  14 

Equilibrium  Volume  Fr actions, (p  ^ »  and  "X  Values  for  Some  Organic  Solvents 


in  Polyethylene. 


Solvent 

32  °F 

<S>2 

70° 

<t>2 

F 

X 

100 

°F 

X 

130 

$  2 

°F 

X 

Water 

0. 0009 

6.  0 

Methanol 

0.  0011 

5.80 

Formic  Acid 

0.  0012 

5.75 

0.0019 

5.28 

0.  0027 

4.  90 

0. 0032 

4.74 

Ethylene  Glycol 

0.  0020 

5.23 

0.0025 

5.00 

0. 0037 

4.  60 

0.0051 

4.28 

Acetic  Anhydride  0.  0025 

5.02 

0.0026 

4.94 

0. 0034 

4.  68 

0.0059 

4.  13 

n-Propanol 

0.0023 

5.  10 

0.  0037 

4.  60 

0. 0062 

4.  08 

0.0099 

3.70 

Diacetone  Ale. 

0.0036 

4.63 

0.  0050 

4.  30 

0. 0072 

3.93 

0.0118 

3.  52 

Nitroethane 

o.  00365 

4.  62 

0.  0052 

4.26 

0. 0067 

4.00 

0.0100 

.3.68 

Glycerine 

0.0050 

4.  30 

0.  0053 

4.24 

0. 0064 

4.08 

0. 0069 

4.02 

n- Butanol 

0. 0046 

4.  37 

0.0066 

4.03 

0.0095 

3.  65 

0.0141 

3.  37 

Aniline 

0. 0064 

4.05 

0. 0090 

3.70 

0.0128 

3.45 

0.  0189 

3.  10 

Butyl  Cellosolve 

0. 0065 

4.  03 

0.  0097 

3.67 

0. 0131 

3.44 

0. 0209 

2.98 

Acetic  Acid 

0. 0066 

4.03 

0.0092 

3.  66 

0.0102 

3.  55 

0.0166 

3.23 

Butyl  Phthalate 

0. 0051 

4.28 

0.0109 

3.  60 

0. 0147 

3.  33 

0.0227 

2.95 

sec -Butanol 

0. 0061 

4.  08 

0.  0115 

3.  55 

0. 0171 

3.20 

0.0283 

2.75 

ter -Butanol 

0. 0085 

3.  76 

0.0130 

3.44 

0. 0189 

3.  11 

0.0310 

2.  67 

Acetone 

0.0108 

3.  62 

0.0144 

3.  35 

0.0183 

3.  13 

0. 0258 

2.83 

n-Octanol 

0.0076 

3.  88 

0.0148 

3.  33 

0. 0189 

3.  11 

0. 0317 

2.65 

Nitrobenzene 

0.0129 

3.45 

0.  0184 

3.  13 

Benz  aldehyde 

0.  0158 

3.27 

0.0214 

3.00 

0.  0288 

2.  74 

0. 0424 

2.41 

ME  Ketone 

0. 0194 

3.08 

0.0252 

2.84 

0.0314 

2.  66 

0. 0447 

2.37 

Ethyl  Acetate 

0. 0216 

2.98 

0.  0272 

2.  77 

0. 0352 

2.57 

0. 0464 

2.33 

Acetaldehyde 

0. 0224 

2.95 

0.0308 

2. 67 

n  -B  utyr  aldehyde 

0.  0246 

2.  87 

0. 0327 

2.63 

0. 0456 

2.  35 

0. 0622 

2.09 

Amyl  Acetate 

0. 0283 

2. 75 

0.0384 

2.48 

0. 0551 

2.  18 

0. 0867 

1.83 

Isopropyl  Amine 

0. 0476 

2.32 

0. 0627 

2.08 

0. 0797 

1.90 

Diethyl  ether 

0. 0882 

1.  88 

0. 0933 

1. 77 

Tetr  adecane 

0. 0910 

1.79 

0. 1030 

1. 70 

0.  1205 

1.  59 

0.  176 

1.34 

Dibutyl  Ether 

0.0800 

1.90 

0.0995 

1. 74* 

0. 1337 

1. 53 

0. 1840 

1.32 

Decane 

0. 0964 

1.75 

0.  117 

1. 62 

0.  157 

1.42 

0.213 

1.23 

n- Heptane 

0.  122 

1.  58 

0.  134 

1. 52 

0.  171 

1. 36 

0.236 

1.  17 

Iso -pentane 

0.  125 

1.  57 

0.  137 

1.  50 

Pentene-2 

0.  143 

1.48 

0.  150 

1.45 

n-Pentane 

0.  133 

1.52 

0.  141 

1.48 

Benzene 

0.  130 

1. 54 

0.  138 

1.  50 

0.  178 

1. 33 

0.237 

1.  17 

p-Chlorotoluene 

0.  128 

1. 55 

0.  146 

1.47 

0.  188 

1. 30 

0.  302 

1. 02 

Chlorobenzene 

0.  134 

1. 53 

0.  147 

1.46 

0. 190 

1.29 

0.  310 

1.01 

Ortho -Xylene 

0.  144 

1.47 

0.  163 

1. 38 

0.228 

1.  18 

0.  319 

0.99 

Par  a -Xylene 

0.  153 

1.43 

0.  167 

1. 37 

0.222 

1.  19 

0.  322 

0.98 

Cyclohexane 

0.  185 

1.  32 

0.208 

1.21 

0.283 

1.06 

0.458 

0.  82 

By  definition,  \  ^  equilibrium  value. 

P 
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17.  Observations  on  the  Master  Solubility  Data. 


The  volume  fraction  §2  recorded  in  Table  14  if  not  equal  to  is  at  least 
equivalent  to  the  solubility  coefficient  S  in  the  relation  P  =  DS.  It  is  not 
surprising,  therefore,  that  the  order  of  increasing  <|>2  in  Table  14  is  es¬ 
sentially  that  already  observed  for  P,  namely:  alcohols,  acids,  nitroderiva- 
tives,  aldehydes,  ketones,  esters,  ethers,  and  hydrocarbons.  Such  inversions 
as  occur  in  the  solubility  order  arise  from  systematic  trends  in  solubility 
for  members  of  any  given  homologous  series,  which  naturally  cause  some 
overlaping  unless  corresponding  members  have  widely  different  solubilities. 

The  P  values  themselves  overlap  less  than  the  solubilities  because  the  dif¬ 
fusion  factor  D  is  also  sensitive  to  polarity,  which  varies  from  functional 
group  to  functional  group;  in  general  D  accentuates  polar  effects  already  dis- 
cernable  in  S. 

The  trend  in  S  within  any  homologous  series  pretty  much  follows  the 
general  solubility  rule  that  "like  dissolves  like'.'.  As  a  polar  functional  group 
of  a  molecule  is  diluted  with  non -polar  hydrocarbon  fragments,  the  molecule 
as  a  whole  approaches  polyethylene  in  all  of  its  properties,  and  the  solu¬ 
bility  increases.  Table  15  illustrates  this  effect.  The  normal  paraffins 
constitute  an  understandable  exception  to  this  rule.  They  are  all  pretty  much 
like  polyethylene  to  start  with  (There  is  from  the  solubility  point  of  view  no 
noticeable  improvement  in  over -all  polarity  in  tetradecane  over  say  pentane.). 
But  it  seems  reasonable  to  expect  that  once  the  polymer  itself  has  set  and 
partially  crystallized- -in  this  partially  crystalline,  partially  amorphous  state 
it  is  believed  that  each  polymer  molecule  may  pass  alternately  through 
several  crystalline  and  amorphous  areas  so  that  in  effect  polyethylene  is  a 
cross-linked  polymer  at  room  temperature,  the  individual  crystallites  serving 
as  tie  points,  -  -reorientation  of  the  polymer  segments  is  restricted.  There¬ 
fore  pentane  is  more  easily  accommodated  than  tetradecane,  for  penetration 
by  the  latter  requires  considerable  rearrangement  on  the  part  of  the  polymer. 
Of  the  solvents  studied,  cyclohexane  seems  to  be  just  about  the  best  in  this 
respect. 
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Table  15 


The  Trend  in  Solubility  Within 
a  Homologous  Series. 


Substance 

Solubility 

70°F 

volume  fraction 

Water 

0.0009 

Methanol 

0. 0011 

n-Propanol 

0. 0037 

n -Butanol 

0.0066 

n-Octanol 

0. 0148 

Formic  Acid 

0. 0019 

Acetic  Acid 

0. 0097 

Nitroethane 

0. 0052 

Nitrobenzene 

0. 0184 

Acetaldehyde 

0. 0308 

n-Butyr  aldehyde 

0.0327 

Acetone 

0. 0144 

ME  Ketone 

0. 0252 

Ethyl  Acetate 

0.0272 

Amyl  Acetate 

0. 0384 

Diethyl  Ether 

0. 0933 

Dibutyl  Ether 

0.0995 

An  interesting  and  significant  feature  of  the  absorption  tests  shows  up 
when  comparing  the  order  of  increasing  solubility  with  the  order  of  in¬ 
creasing  Ep.  Since  a  large  solubility  suggests  a  large  P  whereas  a  large 
Ep  implies  a  small  P,  one  expects  a  large  solubility  to  be  associated  with 
a  correspondingly  small  Ep.  With  respect  to  the  homologous  series  as  a 
whole,  this  is  by  and  large  true.  There  is,  of  course,  a  certain  amount  of 
overlap  to  be  expected  between  the  large  members  of  one  series  and  the 
small  members  of  the  series  above  (on  the  Ep  scale).  Thus  tetradecane  over¬ 
laps  ethyl  acetate,  but  in  general  the  hydrocarbons  with  the  largest  solu¬ 
bilities  have  the  smallest  Ep's  and  the  alcohols  with  the  smallest  solubilities 

have  the  largest  E  's. 

P 

However,  within  any  one  homologous  series  just  the  reverse  tends  to 
be  the  case:  large  solubilities  tend  to  be  associated  with  large  rather  than 

WADC  TR  53-133  Pt  2  39 


small  activation  energies;  the  order  of  increasing  Ep  for  the  straight  chain 
alcohols  is  just  the  order  of  Table  15  and  likewise  for  formic  and  acetic 
acid. 


This  seemingly  anomolous  behavior  is  more  apparent  than  real  when  one 
realizes  that  diffusion  is  much  more  structure  sensitive  than  solubility.  (The 
Ep's  for  ortho-  and  para-xylene  are  15.6  and  13.4  respectively,  although 
their  solubilities  are  nearly  identical.)  An  increase  in  the  length  of  the 
hydrocarbon  chain  attached  to  a  polar  functional  group  almost  invariably  seems 
to  increase  Ep  (the  0.  0348V  term  of  equation  10-2e)  and  just  as  invariably 
increases  the  solubility  in  polyethylene  too.  The  reason  for  the  Ep  increase 
with  size  has  been  discussed.  The  reason  for  the  increase  in  the  solubility 
of  a  polar  molecule  with  increasing  length  of  the  hydrocarbon  chain  is  due  at 
least  in  part  to  the  following  effect. 

A  flexible  hydrocarbon  tail  in  a  polar  molecule  can  swing  back  and 
screen  the  polar  group  from  interactions  with  similar  groups  on  other 
molecules.  In  such  instances  the  molecule  will  appear  less  polar  than  it 
might  otherwise.  The  probability  of  such  internal  shielding  acts  therefore  to 
increase  the  similarity  between  penetrant  and  polyethylene  and  consequently 
enhances  the  solubility. 

It  seems  very  likely  that  the  solubility  data  indicate  the  operation  of 
this  internal  shielding  effect  in  the  series  water,  methyl  alcohol  (no  shielding 
at  all  possible  in  these  two),  n-propyl,  n-butyl,  and  n-octyl  alcohol.  As  one 
should  expect,  the  degree  of  shielding  increases  with  hydrocarbon  length. 

The  molecular  models  also  show  that  acetone  is  an  excellent  example 
of  a  molecule  where  no  internal  shielding  is  possible.  This  helps  explain  its 
position  in  the  tables. 

Shielding  effects  also  help  account  for  the  relative  behaviors  of  the 
three  butyl  alcohol  isomers.  The  effect  is  least  for  n-butanol  with  its  termi¬ 
nal  hydroxyl  group  and  greatest  in  ter -butanol  where  the  polar  hydroxyl 
group  is  embedded  in  a  rather  compact  hydrocarbon  matrix.  These  function¬ 
al  group  structural  features  are  reflected  by  both  their  polarity  factors  (n- 
butanol  :  6.0,  sec-butanol  =  5.4,  ter-butanol  =  4.8)  and  their  solubilities 
(0.0066,  0.0115,  and  0.0130  respectively). 

To  summarize,  then,  as  one  goes  from  say  pentane  to  an  ester  of 
essentially  the  same  shape  and  size,  the  polarity  of  the  molecule  (its  A  AH) 
increases,  the  solubility  decreases  and  Ep  increases.  If  now  one  adds  to 
the  hydrocarbon  portion  of  the  ester,  the  overall  size  of  the  molecule  in¬ 
creases  of  course,  but  its  polarity  decreases,  and  the  solubility,  which  de¬ 
pends  primarily  upon  polarity,  now  begins  to  increase;  E  ,  which  depends 
upon  both  size  and  polarity,  continues  to  increase  although?"* the  two  effects 
work  in  opposition  here  and  seem  nearly  to  cancel  in  the  case  of  the  straight- 
chain  alcohols  whose  Ep  values  all  hover  pretty  closely  around  22.3  kcal/ . 
mole. 
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18.  An  Introductory,  Thermodynamical  Discussion  of  the  Flory-Huggins 
Polymer  Interaction  Parameter  X  * 

The  concept  of  molecular  disorder  of  a  system  and  its  measure,  the 
so-called  entropy  S  has  been  considered.  At  absolute  zero  all  systems  strive 
to  minimize  their  internal  energy  E  (with  a  small  term  PV  added  to  allow 
for  pressure  effects;  i.e.,  actually  it  is  E+PV  -  H,  the  enthalpy,  that  is 
minimized).  The  entropy  of  a  perfectly  ordered  crystal  at  this  temperature 
is  zero.  Observe,  however,  that  although  essentially  all  molecular  motion 
ceases  at  absolute  zero  so  that  there  is  no  thermal  contribution  to  the  over¬ 
all  randomness  or  entropy  of  the  system,  supercooled  liquids  (i.e.,  glasses) 
do  have  a  residual  entropy  at  absolute  zero  relative  to  perfectly  oriented 
crystals  simply  because  of  the  random  orientation  of  their  molecular  con¬ 
stituents. 

We  distinguish,  therefore,  two  types  of  entropy:  namely  entropy  arising 
from  the  random  thermal  vibrations  and  rotations  of  atoms,  molecules,  and 
polymer  segments,  symbol  S^,  and  entropy  arising  from  randomness  in 
position  called  entropy  of  configuration,  symbol  Scf.  The  total  entropy  of 
any  system  is  the  sum  of  these  two. 

Now,  at  any  temperature  T  above  absolute  zero  all  systems  tend  to 
adopt  a  certain  amount  of  disorder:  atoms  begin  to  vibrate,  molecules  begin 
to  rotate,  and  pure  phases  begin  to  dissolve  in  each  other  to  form  random 
mixtures.  This  all  increases  the  entropy  and,  of  course,  the  energy  too, 
for  it  requires  some  of  this  to  do  all  of  these  things.  Clearly  the  energy  is 
no  longer  a  minimum  under  these  conditions,  nor  is  the  entropy  (i.e.,  the 
disorder)  a  maximum,  for  we  could  increase  the  haphazardness  of  the  system 
by  breaking  down  all  molecules  into  atoms,  and  then  these  into  free  electrons, 
protons,  and  neutrons,  etc.  But  at  any  reasonable  temperature  these  later 
possibilities  just  seem  energetically  unlikely.  This  assessment  of  the  situ¬ 
ation  is  reduced  by  thermodynamics  to  the  statement  that  at  constant  T  and 
P  all  systems  tend  to  strike  a  happy  medium  between  maximizing  their 
entropy  (with  a  rise  in  internal  energy)  and  minimizing  their  internal  energy 
(with  a  decrease  in  entropy)  by  actually  minimizing  their  so-called  free 
energy,  F  «  H  -  TS. 

From  this  statement  about  minimizing  the  total  free  energy  of  a  system, 
it  follows  that  the  free  energies  of  the  components  of  a  system  at  equilibrium 
must  be  the  same  throughout  for  otherwise  in  transfering  a  small  amount  of 
one  component  from  one  part  of  the  system  to  another  where  its  free  energy 
is  different  the  total  free  energy  of  the  system  would  change;  if  the  overall 
change  were  negative,  the  process  would  be  a  natural  one  since  it  is  tending 
to  decrease  the  free  energy  and  should  in  time  take  place  spontaneously, 
which  is  to  say  that  the  system  was  not  initially  in  equilibrium.  On  the 
other  hand,  if  the  change  should  turn  out  to  be  positive,  just  the  opposite 
transfer  of  matter  would  lead  to  a  decrease  in  the  total  free  energy  as  before. 

We  can  perhaps  see  from  this  that  the  free  energy  of  a  substance  is 
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essentially  a  measure  of  its  stability.  The  greater  the  free  energy  (large 
H  and/or  small  S),  the  more  likely  it  is  that  the  substance  can  possible  re¬ 
vert  to  some  form  with  a  smaller  free  energy,  i.  e.  ,  the  less  its  stability. 
This  is  sometimes  put  in  another  way:  The  greater  the  free  energy  of  a 
substance,  the  greater  its  escaping  tendency  (as  measured  often  by  its  vapor 
pressure),  activity,  or  so-called  fugacity.  Thermodynamically  these  all  mean 
essentially  the  same  thing. 

For  a  pure  organic  solvent  in  equilibrium  with  a  strip  of  polyethylene, 
this  discussion  informs  us  that  the  escaping  tendency  (say  the  vapor  pressure) 
of  the  solvent  from  the  pure  solvent  phase  must  be  equal  to  the  escaping 
tendency  of  the  solvent  from  the  swollen  polymer  phase. 

Thermodynamically  this  simply  means  that  at  equilibrium 

(18-1)  F^iq  =  F? 

where  as  before  2  stands  for  the  penetrating  species  and  p  for  the  polymer. 
The  little  superscript  bars  merely  indicate  that  it  is  the  actual  contribution 
that  species  2  makes  to  either  the  liquid  or  polymer  phase  that  is  important. 
Sometimes  the  combination  F  is  replaced  by  the  single  symbol ^(mu).  (Also, 
when  using  the  F  notation  the  escaping  tendency  is  moBt  often  referred  to 
as  the  "partial  molal  free  energy"  and,  when  using/* ,  as  the  "chemical  po- 
tential"--all  just  names  and  symbols  for  the  same  thing  however.) 

Writing  out  the  fundamental  equality  above  in  more  detail,  we  get 

(18-2)  flL  -  TSL  =  HP  -  TSP 

2  2  2  2 

which  upon  rearranging  becomes 

(18-3)  SP  -  HL  .  T(SP  -  SL).  or 

2  2  2  2 


(18-4)  ae2  =  tas2 

where  in  general  X  =  (the  X  value  for  the  polymer  phase)  -  (the  X  value 
for  the  liquid  phase). 

These  equations  show  that  from  the  thermodynamic  point  of  view  the 
problem  of  absorption  equilibrium  is  one  of  determining  when  the  heat  of 
solution  of  the  solvent  in  the  polymer,  is  equal  to  the  entropy  of 

solution  (multiplied  by  T),  TAS  . 

w 

According  to  Hildebrand  (Hildebrand  and  Scott,  "The  Solubility  of  Non- 
Electrolytes",  3rd  ed.  ,  pl29),  the  heat  of  solution  is  essentially, 
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(18-5) 


where  the  terms  in  parentheses  raised  to  the  1/2  power  are  the  "solubility 
parameters"  or  "cohesive  energy  density". 

Observe  that  AH^  depends  upon  the  polymer  volume  fraction  <j)p  squared. 
As  this  approaches  zero  (case  of  infinite  swelling,  i.e.,  complete  solution 
of  the  polymer),  so  doesAH2  as  we  might  expect,  for  this  corresponds  to 
mixing  pure  solvent  with  essentially  pure  solvent. 

On  the  other  hand,  for  the  polar,  rather  insoluble  solvents  the  polymer 
volume  fraction,  (1  -Cj)-),  *s  ©8Sentia^Y  unity  even  at  the  point  of  maximum 
swelling  (Table  14).  This  means  that  it  must  be  the  TA§2  term  that  varies 
as  equilibrium  is  being  established  until  finally  equation  (18-4)  is  satisfied. 

At  this  point  it  is  convenient  to  refer  to  our  earlier  discussion  and  de¬ 
compose  AS  2  into  two  parts; 

(18-6)  A§2  =As|h4-  AS§f. 

The  first  term  represents  the  entropy  changes  that  occur  during  mixing 
(solvent  dissolves  in  polymer)  owing  to  alterations  in  the  random  thermal 
motions  (rotations  and  vibrations)  of  the  atoms  and  molecules. 

The  second  term  represents  always  an  increase  in  entropy  that  arises 
upon  mixing  owing  to  the  increased  number  of  configurational  arrangements 
possible  for  the  mixture  over  those  available  to  the  pure  components.  The 
now  classical  statistical  calculation  of  this  configurational  entropy  of  mixing 
for  the  case  of  simple  structureless  systems  dissolving  in  similar  systems 
is  given  in  Appendix  A  to  this  report.  The  problem  becomes  considerably 
more  complex,  however,  when  one  of  the  components  of  the  mixture  is  a 
polymer  composed  of  long  flexible  molecules.  For  one  must  allow  for  the 
fact  in  computing  the  total  number  of  possible  geometrical  arrangements  that 
each  polymer  molecule  may  possibly  be  curled  up,  stretched  out  straight, 
bent  in  a  zig-zag  fashion,  etc.  although,  of  course,  no  two  segments  can  oc¬ 
cupy  the  same  site  (the  so-called  excluded  volume  effect).  These  possibili¬ 
ties  all  contribute  to  an  entropy  of  mixing  that  is  much  greater  than  what 
one  would  obtain  if  the  polymer  molecules  were  treated  as  simple  monomers. 
On  the  other  hand,  the  entropy  of  mixing  is  less  than  what  one  would  calcu¬ 
late  if  each  segment  of  a  polymer  molecule  were  treated  as  an  independent 
species.  Using  the  simplest  possible  model,  Flory  and  Huggins  obtained  the 
following  now  almost  classical  formula  for  AS!** 

(18-7)  AS|*  = '  loge(l-<|>  )+  (For  uncrosslinked, 

^  ^  amorphous  polymers) 

At  this  point  the  theory  stops.  There  is  not  yet  a  reliable  theoretical 
estimate  of  AS^1.  For  non-polar  systems  .thermal  randomness  may  be 
pretty  much  undisturbed  by  mixing,  in  which  case  AS jj1  =  0.  But  only  by 


WADC  TR  53-133  Pt  2 


43 


chance  would  we  expect  AS  to  be  zero  for  polar  solvents  dissolving  in 
polyethylene.  Furthermore,  equation  (18-5)  is  theoretically  valid  only  for 
relatively  non-polar  systems.  In  view  of  all  this  indeterminateness,  the  two 
terms  A  and  TA'Sjj1,  which  appear  together  always  in  the  combination  A  H2 - 
TAS^1,  are  lumped  into  a  single  parameter  "X.  (chi)  multiplied  by  ^  (since  by 
(18-5)  it  is  known  that  in  some  cases  at  least  the  first  term,  AH2,  varies  in 
this  way  with  composition). 

(18-8)  X(RTfp)  =  AH2  -  TA§2h,  or 


X  =  ah2  .  As£h 

RTj)2  TCf 


(18-9) 


(R  :  gas  constant) 


If  we  use  (18-5)  with 
(18-10)  B  =  v°  ^AE^P/v®) 


1/2  -  i^;ap/vo,I/2] 


equation  (18-9)  becomes 
(18-11) 


X--  -  *S*2h 


RT 


1C  is  formally  known  as  the  Flory-Huggins  polymer  interaction  parameter. 


Substituting  now  from  (18-8),  (18-7)  and  (18-6) 
that  at  equilibrium 

(18-12)  AH2-TAS2=  (AH2-TAS2h)  -  TAS  =  (RT<j>^ 


into  (18-4),  we  find 
-f  RT  [loge(l-4p)+^  0 


which  shows  that  the  value  of  the  interaction  parameter  for  swollen  systems 
can  be  obtained  from  solubility  data. 


19.  The  Evaluation  of  the  Polymer  Interaction  Parameter  From 
Equilibrium  Swelling  Measurements. 

From  equation  (18-12)  the  value  of  the  interaction  parameter  at  equili¬ 
brium  is 


(19-1) 


x--  -lojfed-V  -t. 


where  the  volume  fraction^  is  assigned  its  equilibrium  value. 
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Rather  than  substitute  from  Table  16-1  =  1  •)  directly  into 

(19-1)  to  obtain  ~X.  values,  the  function  on  the  r^ht-hand  side  of  equation  (19-1) 
was  plotted  against  the  penetrant  volume  fraction^  as  shown  schematically 
in  Figure  9  and  the  values  of  X  at  equilibrium  read  off  from  this  using  the 
calculated  equilibrium  volume  fractions  as  tabulated  in  Table  14.  The  values 
so  obtained  have  been  included  in  Table  14. 


Fig.  9.  Curve  Used  to  Obtain  "X  at  Equilibrium 


20.  XL  Values  for  Organic  Solvents  and  Polyethylene. 

Although  equilibrium  swelling  measurements  have  been  used  to  obtain 
estimates  of  the  magnitude  of  the  interaction  parameter  (equation  19-1),  it 
is  perhaps  more  natural  to  think  of  X  as  determining  the  degree  to  which 
swelling  may  occur.  As  figure  9  shows,  there  is  a  regular  decrease  in 
the  equilibrium  penetrant  volume  fraction  as  "X  increases. 

This  decrease  in  swelling  with  increase  in  X  is  not  unexpected;  equation 
(18-9)  shows  that  ~K  depends  directly  upon  the  heat  of  solution,  AH2,  of  the 
penetrant  in  the  polymer,  and  swelling  should  decrease  as  AH2  increases; 
since  the  latter  depends  upon  the  polarity  of  the  penetrant,  the  interaction 
parameter  should  show  some  correlation  with  our  polarity  index  AAH  (  Cf. 
Appendix  B). 

Comparison  of  Tables  6  and  14  reveals  that  in  a  general  way  the  order 
of  decreasing  AAH  and  decreasing  X  is  the  same,  namely:  alcohols,  the 
esters  and  ketones,  the  ethers,  and,  finally,  the  hydrocarbons.  Indeed,  Table 
16  demonstrates  that  not  only  is  the  order  of  the  functional  groups  preserved 
in  passing  from  one  parameter  to  the  other,  the  trend  within  any  homologous 
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series  is  preserved  also.  This  is  perhaps  more  than  might  be  expected  at 
first  since  (equation  18-9)  the  interaction  parameter  is  a  free  energy  term 
(actually  the  free  energy  of  mixing  exclusive  of  configurational  effects  di¬ 
vided  by  RT)  whereas  AAH  is  strictly  a  heat  term.  As  we  shall  see  later, 
however,  thgre  is  reason  to  believe  often  that  the  thermal  part  of  the  entropy 
change,  AS2,  varies  with  the  heat  term,  increasing  as  AH2  increases  and 
decreasing  as  decreases.  To  the  extent  that  this  is  so,  X  will  vary  di¬ 

rectly  with  AAH.  Depending  upon  one's  feelings  about  entropy,  one  may  thus 
either  consider  the  correlations  of  Table  16  to  be  "explained",  or  one  may 
prefer  to  consider  the  data  of  the  table  as  strong  supporting  evidence  for  the 
supposition  that  AS varies  with  AH  2  (or  else  is  essentially  constant). 


Table  16 

The  Dependence  of  the  Interaction  Parameter 
Upon  Polarity. 


Substance 

AAH 

kcal/mol. 

-v 

70°F 

Water  (approx.) 

9 

6.00 

Methyl  Alcohol 

7.  1 

5.80 

n-Propyl  Alcohol 

6.3 

4.60 

n-Butyl  Alcohol 

6.0 

4.03 

sec -Butyl  Alcohol 

5.4 

3.  55 

ter -Butyl  Alcohol 

4.8 

3.44 

n-Octyl  Alcohol 

4.  8 

3.  33 

Acetic  Acid 

6.  1 

3.66 

Acetone 

3.0 

3.  35 

Methyl  Ethyl  Ketone 

2.8 

2.84 

Benz  aldehyde 

2.5 

3.00 

n-Butyr  aldehyde 

2.4 

2.63 

Ethyl  Acetate 

2.8 

2.77 

Amyl  Acetate 

2.  5 

2.48 

Isopropyl  Amine 

2.  1 

2.08 

Diethyl  Ether 

1.2 

1.77 

Dibutyl  Ether 

0.  8 

1.74 

n-Decane 

0.0 

1. 62 
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The  following  features  of  Table  16  are  of  special  interest  at  this 
point.  First  of  all, 

a.  It  was  observed  that  ethyl  acetate  and  methyl  ethyl  ketone  have 
similar  P-Factors  at  each  of  the  temperatures  studied;  their  solubilities  in 
polyethylene  are  also  nearly  the  same  all  along  the  line.  Of  the  two,  the 
acetate  has  the  larger  solubility  and  the  smaller  Ep.  Both  of  these  facts  may 
reflect  the  possibility  that  exists  for  greater  screening  in  the  acetate  of  the 
ketone  oxygen  (the  center  of  polarity  in  the  molecule)  by  the  ethyl  group. 


Table  17 

P-Factors  and  Solubilities  of  Ethyl  Acetate  and 
Methyl  Ethyl  Ketone 

_ (From  Tables  V  and  XXIV  of  Part  1) _ 

Ethyl  Acetate  Methyl  Ethyl  Ketone 

Temp  °F  P-Facto^^  -  Solubility ^  P -Factor^  ^  Solubility^ 


32 

1.9 

0.0218 

3.7 

0.0175 

70 

16.  5 

0.0273 

12. 6 

0.0225 

100 

83.3 

0.0353 

60.0 

0.0282 

130 

378. 

0.0466 

326.  6 

-- 

165 

1703. 

•  — 

1401. 0 

_  _ 

(1)  gms/24  hrs/0.001  in.  / 100  in.  ^ 

(2)  gms/gms  polyethylene 

b.  Two  amines  were  studied  in  this  project:  aniline  and 
isopropyl  amine.  The  P-Factors  for  the  former  do  not  fall  on  a  smooth 
curve  on  the  Master  Plot,  so  no  E  value  has  been  calculated  for  it.  And 
only  two  points  are  available  for  the  other,  which  boils  at  93°F.  The  E 
value  calculated  from  these  is  3  kcal/mole  greater  than  that  obtained  usiSg 
equation  (10-2e).  This  all  leaves  the  amine  question  an  open  one  although 
there  seems  to  be  no  compelling  reason  to  expect  them  to  deviate  greatly  from 
the  norm  established  by  the  esters  and  ketones  of  similar  polarity.  In  this 
respect,  it  seems  not  unusual  to  find  the  solubility  (or  X  value)  of  isopropyl 
amine  in  line  with  the  other  compounds  listed  near  it  in  Table  16. 

c.  The  regularity  of  theAAH  -X  correlation  implies  that  the  latter 

can  be  estimated  from  the  former.  Since  one  can  estimate  AAH  values  from 
Dreisbach's  tables  knowing  only  boiling  points,  this  -  X  correlation  means 

that  equilibrium  up-take  values  can  be  predicted  from  readily  available  physi¬ 
cal  data.  This  assumes,  of  course,  that  the  character  of  the  polyethylene  is 
the  same.  Even  if  it  weren't,  the  order  of  Table  16  should  be  unaltered. 
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Observe  that  for  substances  with  uA H  values  between  2  and  3  the  X  value  at 
70°F  is  just  about  numerically  the  same  as  AA.H. 

d.  Finally,  it  is  probably  significant  that  the  hydrogen-bonded 
substances  (the  alcohols  and  acids)  exhibit  Xvalues  noticeably  smaller  than 
one  would  expect  from  their  values.  That  is,  with  regard  to  solubility 

at  least  they  don't  look  as  polar  as  their  heats  of  vaporization  would  suggest. 
Acetic  acid,  which  is  known  to  have  strong  dimerization  tendencies,  is  es¬ 
pecially  anamolous,  being  slightly  less  soluble  in  polyethylene  than  acetone. 
Apparently  association  through  hydrogen  bond  formation  largely  nullifies  the 
inherent  polarity  of  this  substance,  as  seen  upon  vaporization,  so  far  as  the 
comparatively  mild  solution  process  is  concerned.  The  substance  dissolves 
as  the  relatively  non-polar  dimer  rather  than  the  highly  polar  monomer. 

So  much  for  Table  16  and  X  values  at  70°F.  Next  to  be  considered  is 
the  temperature  variation  of  the  interaction  parameter.  Under  favorable 
conditions  such  information  ^ghould  in  principle  (equation  18-9)  yield  infor¬ 
mation  on  both  AH  2  and  AS  ->  . 


2 1 .  The  Temperature  Variation  of  the  Polymer  Interaction  Parameter. 

The  form  of  equation  (18-11)  suggests  a  plot  of  the  interaction  parame¬ 
ter  against  1 /T  at  constant  ^p.  Unfortunately  this  is  not  possible  since  X  is 

only  known  for  the  swollen  polymer  systems  at  equilibrium,  and  the  volume 
fraction  at  equilibrium  naturally  changes  from  one  temperature  to  the  next. 
Nevertheless,  for  the  more  polar  substances  the  absolute  values  of  $2  31,6 
small  and  hence,  also,  their  variations;  i.e.,  for  the  slightly  soluble 
substances  (p  p  is  always  nearly  unity. 

Plotting  X  values  for  the  polar  compounds  against  1/T,  one  finds  that 
the  70,  100,  and  130  degree  points  generally  fall  on  or  near  a  straight  line 
but  that  the  32<^Fpoint  almost  invariably  indicates  an  abnormally  high  solu¬ 
bility  at  this  temperature.  On  the  assumption  that  the  hygroscopic  nature  of 

the  penetrants  was  at  least  in  part  responsible  for  this  anomolous  behavior, 

the  following  data  were  obtained  from  the  slopes  and  intercepts  at  1/T  -  0 
of  the  best  straight  lines  through  the  70,  100,  and  130  degree  points. 
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Table  18 


Approximate  Heats  and  Entropies  of 
Solution  in  Polyethylene. 


Substance 

AH, 

kcal 

A§'h 
e.  u. 

n-Propanol 

4.  8 

7.2 

Dibutyl  Phthalate 

4.  3 

7.4 

Ethylene  Glycol 

4.  1 

4.0 

n-Octanol 

4.0 

6.9 

n -Butanol 

3.8 

4.  8 

sec,  ter-Butanol 

3.9 

6. 1 

Formic  Acid 

3.  3 

0.5 

Acetic  Acid 

2.  6 

1.5 

Acetone 

3.  0 

3.  5 

The  values  of  Table  18  are  probably  accurate  to  only  one  significant 
figure,  if  even  that.  The  relative  values  scaled  to  acetone  are  presumably 
better,  and  it  seems  fairly  certain  that  the  process  of  solution  for  the  car¬ 
boxylic  acids  is  somewhat  unique  considering  the  great  contribution  of  the  car¬ 
boxyl  group  to  the  cohesive  energy  of  organic  liquids  (Table  IV).  In  view  of 
the  inherent  large  polarity  of  the  carboxylic  acids,  the  relatively  small  AR2 
and  (especially)  indicate  that  in  these  cases  the  solute  exists  as  a  dimer 

in  polyethylene. 

Another  way  to  express  this  argument  would  be  to  observe  that  AH 2 
should  be  given  by  »aAH,  and  indeed  this  is,  perhaps  somewhat  fortuitously, 
exactly  true  for  acetone.  But  for  all  the  other  substances  listed  in  Table  18 
AH2<"v  AAH;  the  heat  of  solution  is  less  than  the  polarity  index  would  suggest. 
The  reason  for  this  is  that  our  polarity  index  AAH  has  been  evaluated  from 
data  for  the  monomer  species  whereas  it  appears  that  it  is  the  less  polar 
dimer  species  that  is,  in  fact,  of  interest  here. 

This  same  problem  arose  in  the  calculation  of  Ep  for  highly  associatiated 
substances.  If  the  dimer  is  the  predominant  diffusing  species,  then  we  should 
use  AAHjjjjpgj.  rather  than  in  the  2.4A&.H  term  of  equation  (10-2e). 

Lacking  information  on  AAH^^gy ,  we  were  unable  to  check  this  suggestion. 

But  now  it  seems  that  we  perhaps  can  in  at  least  one  instance. 

The  permeability  data  and  the  temperature  variation  of  the  permeability 
together  with  the  solubility  data  and  the  temperature  variation  of  the  solubili¬ 
ty  data  all  strongly  suggest  that  in  the  case  of  acetic  acid  and  formic  acid  the 
dimer  is  the  principle  solute  species  in  swollen  polyethylene,  even  during  dif¬ 
fusion.  For  these  carboxylic  acids,  then,  it  is  &/\Hdirner  that  we  require  in 
equation  (10-2e).  Moreover,  it  is  just  this  quantity  that  is  reflected  by  the 
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heat  of  solution,  which  we  have.  The  chain  of  reasoning  is  thus  this: 


(21-1) 


*'>Hfor  acetic  acid 
and  polyethylene 


:AAH(^rner  =  AS 2  -  2.6  kcal/.mole 


rather  than  the  value  entered  in  Table  6,  6.  1  =  AAH  .  Using  this 

dimer  value  in  equation  (10-2e)  and  remembering  that  for  a  linear  dimer  the 

size  factor  is  twice  that  for  the  monomer  and  the  shape  factor  is  essentially 

the  same,  we  calculate  for  E  a  value  of  17.1  as  against  the  observed  value 

of  17.8  and  the  monomer  value  of  23.4  kcal/mole.  For  formic  acid  the 

values  are  (Eca*c)  =  20.2,  (Ecalc\  *  15.9,  and  E°bserved -s 

'  p  '  monomer  '  P  'dimer  p 

16.7.  In  each  case  our  improved  calculation  is  low  by  about  0.7  to  0.9  kcal/ 

mole.  This  is  to  be  expected  if  a  small  part  of  the  acid  diffuses  as  the 

monomer . 


There  are  several  factors  that  make  it  difficult  to  extend  this  sort  of 
calculation  to  the  alcohols.  The  strength  of  the  hydrogen  bond  is  less  than  in 
the  carboxylic  acids  and  consequently  the  degree  of  dissociation  upon  diffusion 
will  be  greater.  Also,  the  shape  of  the  dimers  may  vary,  and  there  may  even 
be  trimers,  etc.  in  some  cases. 

When  we  come  now  to  consider  the  solubility  behavior  of  less  polar, 
more  soluble  substances  such  as  the  hydrocarbons,  a  number  of  new  compli¬ 
cating  features  enter  the  picture.  Here  for  really  the  first  time  it  appears 
that  we  must  in  some  way  explicitly  allow  for  the  structure  of  polyethylene 
itself. 


22.  •  The  Solubility  Behavior  of  the  Hydrocarbons  and  the  Structire  of 

Polyethylene. 

Perhaps  the  very  first  thing  one  observes  regarding  the  interaction 
parameter  -  1/T  plots  is  their  marked  departure  from  linearity.  Scarcely 
ever  is  it  possible  to  fit  even  three  of  the  points  with  a  straight  line  whereas 
the  32°,  70°,  100°  and  130°F  points  all  often  lie  on  a  smooth  curve,  the  curve 
being  such  that  the  solubility  seems  to  increase  too  rapidly  with  increasing 
temperature.  This  seemingly  anomolous  hydrocarbon  solubility  behavior  is  at 
first  inspection  most  disconcerting  but  upon  reflection  rather  revealing.  At 
this  writing  the  data  have  not  been  quantitatively  accounted  for.  It  seems 
very  likely,  though,  that  any  satisfactory  explanation  of  the  hydrocarbon 
problem  will  be  mainly  a  matter  of  giving  quantitative  expression  to  the  follow¬ 
ing  observations. 

a.  The  reduction  of  equation  (18-9)  to  (18-11)  depended  upon  the 
fact  that *s  proportional  to<|>p  (equation  18-5).  And  even  if  AH 2  isn't 
given  precisely  by  (18-5)  nor  B  by  (18-10),  the  latter  may  still  be  nearly 
independent  of  the  polymer  volume  fraction  over  a  considerable  range  so  that 
the  form  of  (18-11)  follows  with  B  now  a  to-be -determined  empirical  constant. 
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,  The  second  term  of  (18-11)  is,  however,  probably  not  independent  of 

y  .  AS**1  must,  of  course,  approach  zero  as  the  polymer  becomes  infinitely 

p  c 

dilute,  but  it  need  not  depend  precisely  upon  the  square  of  the  polymer  volume 
fraction  to  do  this,  and  if  the  dependence  upon  is  anything  less  than  this, 
the  whole  term,  AS^/R^  ,  will  increase  as  the  polymer  volume  fraction  de- 
creases  with  increasing  Femperature,  more  will  thus  be  subtracted  from 
the  positive  B /RT  term  at  high  temperatures  than  at  low,  and  curvature  in  the 
Xvs  1/T  plot  will  tend  to  be  as  observed.  In  short,  if  the  entropy  term  is 
not  proportional  to  the  square  of  the  polymer  volume  fraction,  the  interaction 
parameter  will  be  a  function  of  the  composition  in  addition  to  the  temperature. 
This  suposition  is  subject  to  direct  experimental  verification,  (equilibrium 
absorption  studies  at  different  partial  pressures). 

b.  Then  there  is  the  partial  crystallinity  of  the  polymer  to  con¬ 
sider.  As  the  temperature  is  raised,  the  material  becomes  progressively 
less  crystalline.  This  raises  problems  of  a  very  fundamental  sort.  For  sup¬ 
pose,  as  seems  likely,  that  the  penetrants  penetrate  only  the  amorphous 
regions;  that  is,  we  suppose  that  there  are  no  solid  solution  effects.  The 
crystalline  material  remains  pure.  Then  of  course  our  solubility  data  should 
be  reported  as  grams  absorbed  per  gram  of  amorphous  polyethylene,  which 
may  be  less  than  50%  of  the  total.  This  is  the  really  fundamental  quantity. 
Unfortunately  it  is  not  easily  obtained  either  by  theory  or  experiment.  Of 
the  usual  experimental  methods  for  determining  percent  polymer  crystallinity, 
probably  only  the  x-ray  method  would  work  here. 

Nevertheless,  it  is  easy  enough  to  see  in  what  direction  our  data  should 
be  altered.  In  every  case  of  partial  crystallinity  the  "actual"  concentration 
of  the  penetrant  is  greater  than  we  first  thought.  From  figure  9  it  is  ap¬ 
parent  that  this  means  that  all  of  the  X. values  in  Table  14  should  be  lowered. 
(For  the  slightly  soluble  polar  substances  and  50%  crystallinity,  the  correction 
is  about  -0.69  in  X  .  Thus  the  AH2  values  given  in  Table  18  will  be  unaf¬ 
fected,  the  full  brunt  of  the  correction  falling  onAS^1.)  We  have  then  that 
the  interaction  parameter  is  essentially  infinite  for  crystalline  polyethylene 
(no  swelling)  and  less  than  the  figures  given  in  Table  14  for  amorphous  poly¬ 
ethylene  (limited  swelling).  (Incidentally,  figure  9  shows  that  swelling  is 
infinite  at  X  equal  to  0.5.) 

Now  we  come  to  an  important  point.  Of  all  the  values  listed  for  a  given 
solvent  in  Table  14,  the  32°  values  will  have  to  be  corrected  the  most  and 
the  130°  values  (less  crystallinity)  the  least.  This  may  help  explain  the  curva¬ 
ture  in  the  apparent  solubility  curves;  as  the  temperature  is  raised,  the  solu¬ 
bility  does  not  actually  increase  as  rapidly  as  the  absorption  tests  suggest  for 
the  effective  weight  of  the  polymer  (i.e.,  amount  of  amorphous  material)  is 
increasing.  But  if  this  is  true  of  the  hydrocarbons,  why  isn't  it  also  true  of 
the  other  solvents  ? 

Well,  to  some  extent  it  is,  and  for  this  reason  the  data  of  Table  18  are 
only  approximate  as  indicated.  But  for  the  slightly  soluble  substances  the 
curvature  in  the  Xvs.  1/T  plots  was  much  less  than  it  was  for  the  less  polar 
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substance.  This  is  perhaps  not  altogether  surprising.  Polyethylene  com¬ 
pletes  its  melting  somewhere  around  240°F.  At  130°  F  there  probably  isn't 
too  much  melting  going  on  in  the  pure  polymer,  and  as  the  volume  fraction 
figures  in  Table  14  show,  polymer  swollen  to  its  maximum  with  say  butyl 
alcohol  or  even  acetone  is  still  nearly  pure  polymer  at  all  the  temperatures 
studied.  But  this  is  far  from  the  case  for  the  more  soluble  hydrocarbons, 
especially  at  130°F.  At  such  high  penetrant  concentrations  the  properties  of 
the  swollen  polymer  are  likely  to  be  significantly  different  from  the  properties 
of  the  pure  polymer  at  the  same  temperature. 

c.  In  particular,  the  crystalline  polymer  -  amorphous  polymer 
equilibrium  will  be  disturbed  toward  the  amorphous  polymer  side.  Another 
way  of  putting  this  is  to  say  that  the  melting  points  of  the  crystalline  regions 
are  lowered  as  one  dissolves  "anti -freeze"  in  the  amorphous  regions. 

In  this  respect  ethylene  glycol  is  a  poor  "anti-freeze"  owing  to  its  very 
limited  solubility  in  amorphous  polyethylene.  On  the  other  hand  the  hydro¬ 
carbons  with  their  inherently  large  solubility  in  amorphous  polyethylene  are 
"good"  anti-freezes  and  may  well  cause  a  significant  amount  of  melting  of 
the  crystalline  regions  at  I30°F. 

With  this  discussion  of  polymer  crystallinity,  freezing  point  lowering 
effects,  and  AS  2  behind  us,  we  are  in  a  position  to  introduce  an  important 
although  heretofore  neglected  feature  of  swelling  equilibrium  in  polyethylene. 

d.  It  may  be  introduced  this  way:  Why  doesn't  polyethylene  dis¬ 
solve  completely  in  a  hydrocarbon  such  as  decane  ?  After  all,  decane  is  a 
npn-polar  substance  (a&H  s  0).  Therefore  its  heat  of  solution,  AH^,  in 
polyethylene  should  be  small  (zero  according  to  our  earlier  discussions  of 
this).  This  implies  (equation  18-9)  a  small  value  for 'X.  and  therefore  an  infi¬ 
nite  amount  of  swelling  (Figure  9.  Polymers  dissolve  in  any  solvent  for 
which  ~X  <  0.5).  Yet  actually  the  swelling  in  decane  is  limited.  This  almost 
forces  one  to  conclude  that  polyethylene  is  in  effect  cross-linked.  (Figure  9 
applies  to  only  non-cross -linked  polymers.  For  cross-linked  polymers  the 
swelling  is  always  limited. )  This  notion  has  already  been  introduced  in 
Section  17.  It  was  pointed  out  there  that  each  polymer  molecule  may  pass 
alternately  through  several  amorphous  and  crystalline  regions.  Thus  no 
matter  how  much  the  amorphous  regions  may  swell  the  polymer  as  a  whole 
remains  rigid  so  long  as  there  are  any  appreciable  number  of  crystalline  tie 
points  present. 

It  may  at  first  seem  strange  that  cross-linking  should  affect  the  inter¬ 
action  parameter  primarily  through  the  entropy  term  of  equation  (18-9)  rather 
than  the  energy  term,  AH But  Flory  and  Rehner  have  shown  that  this  is 
the  case.  The  physical  picture  is  simply  this:  as  a  cross-linked  polymer 
is  swollen  by  a  penetrant  the  segments  of  the  polymer  molecules  in  the 
swollen  regions  are  stretched  out  into  unnatural  elongated  shapes.  This  de¬ 
creases  their  configurational  entropy.  The  theory  suggests  a  simple  cor¬ 
rection  of  (-Rv2  ^  p  ^/vc)  to  (18-7),  the  fundamental  equation  for  non-cr.oas- 


WADC  TR  53-133  Pt  2 


52 


linked  polymers.  Here  as  before  v°  is  the  molar  volume  of  the  pure  pene¬ 
trant;  vc  is  the  so-called  molar  volume  of  the  polymer  between  cross-links. 
For  a  polymer  in  which  there  is  no  cross-linking  vc  is  infinite  and  the  cor¬ 
rection  term  above  vanishes. 


This  all  means  that  for  cross-linked  polyethylene  (i.e.,  any  crystalline 
polyethylene  of  average  molecular  weight)  the  term  AS^1  we  introduced  as 
the  purely  thermal  part  of  the  entropy  of  mixing  does  in  fact  contain  a  con¬ 
figurational  part,  namely  the  entropy  of  elastic  stretching  of  the  polymer. 
Using  the  Flory-Rehner  expression  for  this,  we  write 

(22-1)  Asf'AS^'  -  Rv2°^/3/vc 


where  now,  insofar  as  the  theory  has  properly  assessed  .configurational  ef¬ 
fects,  AS^'  is  a  purely  thermal  entropy  term. 


Equation  (22-1)  provides  a 
non- vanishing  X.  values.  Even 


a  quantitative  explanation  of  the  hydrocarbon 
if  AH  2  -  AS^'  =  0,  one  still  has  that 


(22-2) 


* 


♦ 


2 

573- 


Equation  (22-2)  suggests  several  interesting  "experiments"  none  of 
which  has  been  tried. 


For  example,  what  values  do  the  data  on  the  various  hydrocarbons  give 
for  vc,  and  how  does  v  vary  with  the  temperature?  Since  v  is  a  parameter 
characteristic  of  the  pofymer  itself,  its  value  should  be  the  same  in  all 
penetrant-polymer  systems  of  the  same  degree  of  crystallinity.  And  reliable 
information  on  the  increase  in  vc  with  temperature  would  tell  us  just  how 
melting  proceeds  in  these  systems.  The  difficulty  of  course  is  that  the 
melting  information  is  required  in  order  to  obtain  vc  in  the  first  place. 

Incidentally,  equation  (22-2)  predicts  the  order  of  "X values  for  pentane, 
decane,  and  tetradecane  actually  observed.  The  smaller  the  molar  volume 
of  the  penetrant,  the  smaller  "X  and  <|>p  must  be  if  vQ  is  to  be  constant. 

This  essentially  completes  our  discussion  of  the  absorption  tests.  We 
conclude  with  a  few  remarks  about  the  polarity  index  AAH  and  its  relation  to 
the  cohesive  energy  density,  (C.E.D.  5  AE^P  )  Also  the  thinkness  measure¬ 
ments  will  be  briefly  reviewed.  Then  we  fake  up  the  last  part  of  the  P- 
Factor  problem,  the  PQ  factor  of  equation  (l-2a). 

23.  The  Polarity  Index  AAH  and  its  Relation  to  the  Cohesive  Energy  Density. 

According  to  equation  (18-5)  the  heat  of  solution  of  a  substance  (com¬ 
ponent  2)  in  polyethylene  (component  p)  is  proportional  to  the  square  of  the 
difference  between  the  square  roots  of  their  respective  cohesive  energy  densi- 
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ties,  (the  molar  heat  of  vaporization  divided  by  the  molar  volume).  By 
slight  rearrangement,  we  obtain 

(23-1)  AH2  =4*  (<AE™V/2  -  /v°  )1/2]  2 

which  is  to  be  compared  to  our  A\H, 

(23-2)  =AE^ap  -  AEvap  (v?  /vo  ) 

2  p  2  p 

where  the  second  term  of  (23-2)  clearly  corresponds  to  the  second  term  of 
equation  (5-3):  "the  heat  of  vaporization  of  a  hydrocarbon  of  similar  size." 

Although  the  two  equations  look  quite  alike  when  is  nearly  one  (the 
case  for  which  AAH  was  introduced),  they  yield  rather  Afferent  values.  For 
example,  if  the  two  terms  in  (23-2)  are  16  and  9  respectively,  r  7 

whereas  by  (23-1)  Ai^  *  L 


24.  Swelling  Index. 

These  data  are  reported  as  the  thickness  of  the  swollen  polyethylene 
strips,  L,,  minus  their  original  thickness,  L°,  divided  by  L°-  Let  us  call 
this  quantity  \ ' 

(24-1)  \=  (L-L°)/L°, 

and  let  Lj,  L>2?  and  Lj  be  the  final  dimensions  (those  of  the  swollen  strip) 
and  Lj,  L£»  and  the  original,  dimensions  (those  of  the  unswollen  poly¬ 
ethylene  itself)  of  a  strip.  Then  in  general 

(24-2)  L-  *  Lf(l-^i)  i  *  1,2,3. 

If  the  strips  are  isotropic,  Xj  *  ^-2  =  ^3  *  justAsay.  This  very  possibly  is 
not  entirely  true  in  the  present  instance,  but  we  shall  have  to  assume  that  it 
is.  For  then  we  can  write  that  the  volume  V  of  the  swollen  polymer  is 

V  =  L!L2L3  .  L°L2Lf(l+X)3 

where,  of  course,  L^L§L§  *  V°  ,  the  volume  of  the  pure  polyethylene  former 
ly  symbolized  as  v°).  Expanding  and  rearranging,  we  find  that 

(24-3)  (V-Y°)/V°  *  3A(1+  ^  +  \2/3). 

This  expression  is  to  be  compared  with  our  assumption  earlier  (Section 
16)  that  the  volume  V  of  the  swollen  polymer  is  just  the  volume  V°  of  the 
pure  polymer  plus  the  volume  m  2.1^2  ° *  pure  penetrant.  (Here  m2* 

weight  of  penetrant  =  increase  in  weight  of  the  strip).  The  volume  additivity 
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assumption  implies  therefore  that 


(24  -4)  (V -V°)/V°  •  (m2  /D2)  /(mp  /Dp)  -  Q(Dp  /D2) 

where  Q  *  m2/mp  ."grams  absorbed  per  gram  of  polyethylene." 

The  two  expressions,  3W  X  *  X2 /3)  and  Q(Dp/D2),  have  been  com¬ 
pared  at  32°  70°,  100°  and  130°F  for  all  substances  as  soluble  as  ethyl 

acetate ;  for  substances  less  soluble  than  ethyl  acetate  the  swelling  was  really 
almost  too  small  to  be  determined  accurately  by  direct  linear  measurement. 
In  fact,  the  weight  measurements  generally  seem  to  be  approximately  three 
times  as  precise  (as  judged  by  their  reproducibility)  as  the  linear  measure¬ 
ments.  Table  19  lists  some  representative  values  for  these  two  expressions. 
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Table  19 


A  Comparison  of  the  Absorption  Tests  and  the  Thickness 
Data.  AV/V°  =  3>j(l-+  X  +  X2/3)- 


Substance 

32°  F 

70°F 

100°F 

130°F 

Ethyl  Acetate 

QPp/Efe) 

0.0221 

0.0280 

0.0365 

0.0486 

AV/V° 

0.01 

0.014 

.015-. 034 

.020-.  048 

Amyl  Acetate 

Q(Dp/r^) 

AV/Vo 

0.0291 

0.0400 

0.0582 

0.0950 

.015-. 024 

.024-. 039 

.026-.  043 

.068-.  093 

Dibutyl  Ether 

QiVp  /d2) 

0.0869 

0. 1108 

0. 1556 

0.225 

AV/V° 

.048-.  058 

.089-. 102 

. 128-. 139 

. 146-. 148 

n-Butyr  aldehyde 

Q(Dp/D2) 

0.0252 

0.0339 

0.0477 

0.0663 

AV/VO 

.015-.  024 

.024-.  034 

.014-.  024 

0.038 

n -Pentane 

Quycfc) 

AV/V° 

0.  154 

0.  153 

0.  164 
.169-.  181 

n -Heptane 

Q(DP/D2( 

AV/V° 

0.  138 

0. 155 

0.207 

0.309 

0.  141 

. 141-. 160 

. 182-. 193 

0.298 

n-Decane 

Q(Dp/D2) 

0.  107 

0.  133 

0.  186 

0.271 

AV/V° 

.098-.  085 

.091-. 098 

. 147-. 151 

0.232 

n-Tetradecane 

/d2) 

0.  100 

0.  116 

0.  137 

0.212 

AV/V° 

.068-.  088 

. 078 -. 086 

. 131- . 149 

. 171-.206 

Benzene 

Q(I^  /d2) 

0.  150 

0.  160 

0.217 

0.  312 

AV/V° 

.  122 -.  138 

. 152-. 172 

.202-. 211 

.244-. 33 

o -Xylene 

Q(4>  /Efe ) 
av/v° 

0.168 

0.  195 

0.295 

0.586 

. 153-. 171 

.  1 84-.  2 12 

.270-. 292 

0.50 

p-Xylene 

Q(Dp/D2) 

AV/VO 

0.  180 

0.200 

0.284 

0.473 

. 145-. 188 

.  162-.  193 

.242-. 278 

.402-. 422 

p  -  Chlor  otolue  ne 

Q(4>  /d2) 

A V/V° 

0.  148 

0.  171 

0.231 

0.433 

.  128-.  148 

. 169-. 199 

. 187-.240 

.  394-.  413 

Cyclohexane 

Q(Efe  /d2) 

AV/V° 

0.228 

.  187-. 194 

0.262 
.266-. 331 

0.394 

0.394 

0.847 
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It  seems  that  about  all  one  can  safely  conclude  from  these  calculations 
is  that  the  actual  swelling  as  determined  by  the  thickness  measurements  is 
generally  equal  to  or  less  than  the  swelling  that  one  calculates  from  the 
additivity  assumption.  That  is,  there  may  be  a  certain  amount  of  filling  in 
of  voids  in  the  polyethylene  by  the  penetrants.  If  the  voids  are  not  equilibri¬ 
um  thermal  ones,  solution  by  this  mechanism  will  increase  the  mass  of  the 
polymer  phase  but  not  its  volume. 

In  principle  measurements  such  as  these  might  under  favorable  con¬ 
ditions  yield  valuable  information  on  polymer -penetrant  interactions.  Unfortu¬ 
nately,  from  this  point  of  view,  the  swelling  of  polyethylene  is  often  scarcely 
detectable.  On  the  other  hand,  in  cases  where  the  swelling  is  easily  measur¬ 
able  the  problem  is  complicated  by  the  "anti-freeze"  effect.  So  we  leave 
these  measurements  now  and  turn  finally  to  the  problem  of  PQ. 


25.  An  Empirical  Treatment  of  the  P0  Data. 

Table  3  reveals  that  there  is  a  general  decrease  in  PQ  with  Ep. 

This  suggests  a  plot  of  the  logarithm  of  P0  (all  logarithms  from  here  on 
are  to  the  base  10)  against  Ep  (in  kcai/mole)  as  shown  in  Figure  10.  Not 
all  the  data  of  T able  3  appear  in  this  figure.  Only  those  homologous  series 
are  represented  for  which  there  are  data  on  at  least  three  members,  namely 
the  aliphatic  hydrocarbons,  the  aromatic  hydrocarbons,  the  aliphatic  acids, 
and  the  aliphatic  alcohols. 

The  extent  to  which  the  nineteen  compounds  plotted  in  Figure  10  group 
themselves  according  to  homologous  series  membership  is  truly  gratifying. 
Furthermore,  the  following  compounds  could  have  been  added  and  still  the 
overall  picture  would  have  been  pretty  consistent. 

Cyclohexane  (with  the  aromatic  hydrocarbons), 

Phenol  (with  the  aliphatic  alcohols),  and 

Dibutyl  ether.  Ethyl  Acetate,  ME  Ketone,  Acetone, 

Amyl  Acetate,  Butyr  aldehyde ,  and  Benz  aldehyde 
(all  with  the  aliphatic  hydrocarbons,  although 
several  of  these  faty  on  a  nice  straight  line  of 
their  own) . 

Such  a  classification  shows  that  the  grouping  in  Figure  10  is  to  some 
extent  determined  by  shape  and,  perhaps  to  a  larger  extent,  by  polarity.  Thus 
non-polar  but  round -shaped  cyclohexane  falls  in  with  the  also  relatively  non¬ 
polar  aromatic  ring  compounds,  and  the  straight  chain  ethers,  esters  and 
ketones  are  grouped  somewhat  with  the  straight  chain  hydrocarbons.  But 
highly  polar  phenol  clearly  belongs  with  the  straight  chain  alcohols,  and  the 
straight  chain  acids  fail  in  a  class  by  themselves  rather  than  with  the  paraf¬ 
fins. 
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a  Aliphatic  VW^rocarbi 

1  -PtWTAHE 

2. n  -HePTANe 

3. n  -  Decoke 

4. n  -  Tetrads  c a* 


+  Aromatic  Hydro cARt 

1.  Beuzeme 

2.  p  -  XyuCHE 

3.  Chloho  eex'z  t 

4.o  -XvutNe 

5.  p  -  CHIOROTOIAH 

x  Auphatic  Acids 

1.  Formic  Acto 

2.  Acetic  Acii 

3.  h-  Bvjtyric  Ac 


©  Auphatic  Alcoho 

4.  Methyu  Aim 

2.  Acuyv.  Awcoi 

3.  h-PRQPYv.  Atco 

4.  ft-  Bvtyl  AlcQi 
S  h~  Octyl  Aico; 
4sec-Bu-nru  Aufcc 
7.  ter-  Butyl.  Auc< 


26.  The  logP^  vs.  Ep  Plot  and  its  Uses. 

Figure  10  in  conjunction  with  our  earlier  discussions  may  in  some 
instances  be  as  or  even  more  useful  than  the  Master  P -Factor  Plot  itself. 

Its  utility  will  be  demonstrated  by  two  .rather  different  examples. 

i  £■  - 

(1)  In  conjunction  with  equa itQf\  Q0~).Q  FitfitTz-  10  can  be  used  to  ob¬ 
tain  absolute  values  for  P -Factors.  It  fJT  tyji/  /}£gessary  to  have  the  usual 
vital  statistics  on  the  compound  in  question  (i.e.,  density,  molecular  weight, 
model  length,  and  heat  of  vaporization)  and  information  regarding  its  family 
ties  in  Figure  10.  For  example,  knowing  that  amyl  acetate  is  more  or  less 
related  to  the  straight  chain  hydrocarbons,  we  find,  using  the  value  19.5  for 
Ep  (Table  13),  that  logPQ  *  15.2  (compared  to  the  observed  value,  Table  3, 
of  15.32).  From  these  values  for  Ep  and  logPQ  the  P-Factor  for  amyl 
acetate  at  70°F  is  calculated  to  be  5.  q.  The  observed  value  was  8.7  (Table 
1). 

It  might  perhaps  be  worth  mentioning  that  a  certain  amount  of  interpo¬ 
lation  is  possible  in  Figure  10  just  as  it  was  in  the  Master  Plot.  Dipentene, 
for  example,  probably  does  not  belong  either  with  the  aromatics  or  the  ali¬ 
phatic  hydrocarbons,  but  rather  some  place  in  between  (see  structural  formula 
in  Section  12).  Just  where  in  between  is  determined  by  its  calculated  Ep 
value,  15.9,  and  guesstimate,  say  13.9  to  14.0  for  logPQ. 

(2)  It  may  be,  though,  that  for  one  reason  or  another  a  trustworthy 
value  for  Ep  is  not  available:  Possibly  the  heat  of  vaporization  is  in  doubt. 

In  this  case  a  single  P-Factor  determination,  say  at  room  temperature,  will 
suffice  to  determine  the  permeability  behavior  at  all  other  temperatures 
provided  the  compound  can  be  assigned  to  one  of  the  families  in  the  logPQ  vs. 
Ep  plot.  We  will  illustrate  the  application  of  Figure  10  to  such  a  problem 
with  allyl  alcohol. 

Let  us  suppose  that  all  we  know  about  allyl  alcohol  is  that  it  is  an  alco¬ 
hol  and  that  at  room  temperature  (70°F)  it  has  a  P-Factor  of  0.65.  What 
is  its  P-Factor  at  165°F? 

If  we  only  knew  E_,  we  could  use  the  70°  P-Factor  and  equation  (l-2a) 
to  find  PQ,  and  then  with  this  and  Ep,  P  could  be  calculated  at  165°F.  This 
sort  of  calculation  has  been  mentioned  before.  Now  this  time,  however,  we 
are  supposing  that  we  do  not  know  E  .  But  we  do  know  that  whatever  Ep 
and  logPQ  may  be  for  allyl  alcohol,  they  must  be  such  that  the  point  for  it 
falls  on  the  alcohol  curve  in  Figure  10.  Of  course  there  are  many  sets  of 
values  for  the  pair  that  satisfy  this  condition,  but  only  one  set  that  will,  when 
used  in  equation  (l-2a),  give  the  right  value  for  P  at  70°F.  So  in  principle 
the  problem  is  solved.  In  practice  it  is  convenient  to  know  the  equation  for 
the  logPQ  vs.  Ep  curve.  From  Figure  10  we  find  that  the  parameters  m 
and  b  in  the  equation 

(26-1)  logPQ  *  mEp  ♦  b 
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have  the  following  values. 


Table  20 

The  Values  of  the  Parameters  in  the  Equation: 


logPQ  =  mEp4  b. 


F  amily 

m 

b 

Aliphatic  Hydcb . 

0.483 

5.  86 

Aromatic  " 

0.  661 

3.74 

Alcohol 

0.  664 

1.40 

Aliphatic  Acid 

1. 28 

-9.  00 

Essentially  what  we  have  then  are  two  unknowns, 
two  relations  between  them: 


Ep  and  logp0> 


and 


logP  ■  -Ep/2.3RT  -t-  logPo,  and 

logP  =  mE„+  b. 

0  o  p 

Adding  and  solving  for  E  ,  we  obtain  in  general 

P 

(26-2)  E  =  (b  -  logP) /( 1/2. 3RT  -  m) 

P 

=  (b  -  logP) /(2 18. 6/T  -  m) 
and  in  particular  for  allyl  alcohol  at  70°F, 

Ep  •  (1.40  -  logO.  65)  /(2 18.  6/294.2  -  0.664) 

s  20.  1  kcal/mole 

compared  to  the  observed  value  of  20.2  kcal/mole.  From  equation  (26-1)  (or 
Figure  10)  it  then  follows  that 

logPQ  *  14.75 

for  allyl  alcohol  (the  observed  value  is  14.85).  Substituting  in  the  logP 
equation  with  T  a  347°K,  we  find  finally  a  value  of  123  for  the  P -Factor  of 
allyl  alcohol  at  165*^.  The  observed  value  is  138.  This  is  remarkable  a- 
greement  but  then  so  is  the  linearity  of  the  curves  in  Figure  10. 

One  of  the  advantages  of  this  method  is  that  if  the  Ep  value  used  is  too 
large,  logPQ  will  also  be  too  large,  and  the  two  errors  will  tend  to  cancel. 
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27.  A  Qualitative  Theoretical  Discussion  of  the  Factor. 

As  we  mentioned  previously,  just  as  the  P-Factors  are  of  the  form 

„  „  -E  /RT 

P  -  PQe  p 

so  too  the  solubility  and  diffusion  coefficients  in  the  relation  P  =  DS  are  of 
the  form 

D  =  D_e"Ed/RT  and 


S  =  S  e-Es/RT 
o 


where 


Ep  * 


¥  Ee  and  P„  =  DR5, 


Actually  the  absorption  tests  data  reveal  that  the  solubilities  of  organic 
liquids  in  polyethylene  are  not  described  as  easily  as  the  simple  exponential 
form  above  would  suggest.  But  for  the  slightly  soluble  substances  the  solu¬ 
bility  behavior  approximates  this  where  Es  is  essentially  AH  2  and  is  small 
compared  to  and  where  SQ  is  essentially  e^S^  / 

In  the  transition  state  theory  (Glasstone,  Laidler  and  Eyring,  "The 
Theory  of  Rate  Processes",  p524)  the  diffusion  constant  receives  a  similar 
treatment.  By  considering  the  equilibrium  between  activated  penetrant  mole¬ 
cules  (molecules  just  in  the  act  of  diffusing)  and  normal  penetrant  molecules 
(ones  at  rest  in  the  polymer),  one  finds  that  E<j  is  the  so-called  "energy  of 
activation  for  diffusion"  and  Dq  is  essentially  e~&S*/R  where  AS*  is  the  entro¬ 
py  change  that  occurs  within  the  system  when  a  normal  penetrant  molecule 
is  "activated.  "  It  is  called  the  "entropy  of  activation". 

Combining  these  statements,  we  have  the  following  equivalence. 

logP0  *  log^o  +  loSs0 

-th' 

OJ  AS*  +&S2  . 

But  we  also  have  the  relation  (equation  26-1) 

logP  *  mE  +  b. 

6  o  p 

It  is  worth  noting  that  in  the  language  of  the  transition  state  theory  E  would 
be  written  AE*.  This  suggests  that  ^ 


mE  ^viAS*  and  hence  also  that 
P 

,  -th' 
b  ro  £  S2  . 
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All  we  need  to  do  is  justify  one  of  these  and  the  other  one  follows.  We 
shall  attempt  to  argue  the  first  one. 

It  says  that  the  entropy  of  activation  is  roughly  proportional  to  the 
energy  of  activation.  That  this  might  indeed  be  so  cam  perhaps  be  seen  from 
these  considerations.  There  are  many  many  ways  a  normal  penetrant  mole¬ 
cule  might  get  activated.  Of  all  of  these,  the  most  probable  one  is  the  one 
with  the  largest  equilibrium  constant  K  ■  (activated  molecules) /{normal  mole¬ 
cules).  (The  polymer  immediately  surrounding  a  molecule  in  the  act  of  dif¬ 
fusion  constitutes  part  of  the  so-called  "activated  complex".  Likewise  the 
polymer  immediately  surrounding  a  molecule  at  rest  constitutes  part  of  what 
we  might  call  a  "normal  complex".  In  the  equilibrium  constant  expression 
above  "activated  complex"  and  "normal  complex"  have  simply  been  written 
'activated  molecules'  and  'normal  molecules'  respectively.) 

Now  from  classical  thermodynamics 

RT logeK(normal  molecules  =  activated  molecules)  =  -AF°* 


where  AF°*  is  the  standard  free  energy  change  for  the  reaction  indicated  a- 
bove  of  normal  molecules  going  to  activated  ones;  i.e.,  it  is  the  so-called 
"free  energy  of  activation".  Clearly  K  will  be  greatest  when  AF°*  is  as 
small  as  possible.  Now,  again  from  classical  thermodynamics, 

AF°*  -  AH0*  -  TAS°*- 

(The  little  superscript  zeros  on  these  terms  merely  indicate  that  configur¬ 
ational  effects  have  been  sorted  out.  For  most  purposes  it  really  isn't  neces¬ 
sary  to  show  this  for  the  AH  term,  but  as  we  have  seen  configuration  can 
have  a  large  effect  upon  the  entropy  and  hence,  as  the  equation  shows,  upon 
the  free  energy  too.) 

In  this  equation  the  AF°*  term,  for  example,  represents  the  difference 
between  the  free  energy  of  an  activated  molecule  and  a  normal  one.  Since 
we  start  with  the  latter,  nothing  can  be  done  in  the  way  of  adjusting  its  free 
energy  in  order  to  minimize  AF0*.  This  leaves  only  the  free  energy  of  the 
activated  molecules  to  examine.  If  we  let  Factivated  moieCule  5  F*’  etc- 

this  means  that  K  will  be  as  large  as  possible  when 

F*  =  H*  -  TS* 


is  as  small  as  possible. 

At  this  point  we  encounter  a  fundamental  aspect  of  nature.  F*  is  not 
minimized  by  making  H*  as  small  as  possible.  Rather,  it  is  minimized  by 
striking  some  sort  of  compromise  between  a  not  too  large  H*  and  a  not  too 
small  S*.  This  argument  may  sound  a  bit  familiar,  for  we  have  used  it  be¬ 
fore --several  times.  K*  could  conceivably  be  large  and  S*  small,  but  this 
is  (thermodynamically)  unlikely,  and  the  converse  can  never  be  true.  So 
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there  are  good  sound  reasons  why  S*  should  follow  H*,  which  is  essentially 
equivalent  to  Ejj.  And  since  E^  is  the  major  part  of  Ep,  we  have  established 
what  we  wished  to  establish. 

This  correlation  of 


mE  &S* 
P 


AS 


activation 


-th' 

b— HS2 


as 


solution 


and 


helps  explain  several  interesting  features  of  Table  20,  namely  the  decrease 
in  b  and  the  increase  in  m  in  going  from  the  aliphatic  hydrocarbons  to  the 
aliphatic  acids.  It  seems  reasonable  to  expect  that  the  entropy  of  solution 
will  be  small  whenever  there  is  a  tendency  toward  dimerization.  This  ac¬ 
counts  for  the  trend  in  b.  The  trend  in  m,  although  not  at  all  surprising,  is 
not  rationalized  so  readily.  A  small  entropy  of  solution  of  course  leaves  a 
lot  to  be  gained  during  activation  if  the  dimers,  or  even  just  a  part  of  them, 
are  split  during  diffusion.  But  this  could  come  about  without  any  increase  in 
the  m  factor  of  the  product  mE„  since  Ep  will  naturally  be  larger  than  it 
might  otherwise  be  if  there  is  dimer  dissociation  during  diffusion. 


Well,  in  any  event  there  are  trends,  and  a  plot  of  m  against  b  is  best 
fitted  by  a  straight  line  connecting '  the  two  extreme  points  (the  aliphatic  hydro¬ 
carbon  and  aliphatic  acid  points)  and  passing  between  the  two  intermediate  ones 
which  don't  lie  too  far  from  it. 


28.  Summary. 

Neither  polyethylene  nor  the  organic  penetrants  of  this  project  should 
be  considered,  even  by  themselves,  as  simple  substances,  for  there  is  not 
yet  a  satisfactory  theoretical  explanation  for  even  the  commonest  physical 
properties  of  this  polymer,  nor  is  the  liquid  state  entirely  understood  even 
for  substances  as  simple  as  liquid  argon.  It  is  not  surprising,  therefore,  that 
the  interaction  of  polyethylene  with  an  organic  penetrant  turns  out  to  be  a  rela¬ 
tively  complex  problem  in  statistical  thermodynamics  and  diffusion  theory.  Now 
no  such  problem  has  ever  been  worked  out  exactly  even  in  the  very  simplest 
possible  cases.  This  report  shows,  nevertheless,  that  it  is  possible  to  draw 
a  number  of  interesting  conclusions  from  the  data  supplied  by  Plax  on  the 
permeability  of  polyethylene  to  organic  solvents.  These  conclusions  are  sum¬ 
marized  below.  Their  theoretical  justification  and  the  exact  conditions  under 
which  they  may  be  expected  to  be  valid,  were  established  in  the  earlier 
sections  of  this  report. 

1.  An  inductive  picture  of  the  permeability  of  polyethylene  to  organic 
solvents  is  best  based  upon  the  relation:  P  =  DS. 

2.  Of  the  two,  S  may  be  taken  to  be  much  more  important  than  D. 
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3.  Very  often  the  order  of  S  from  swelling  measurements  is  also 
precisely  the  order  of  the  observed  P -Factors. 

4.  The  solubility  coefficient  S  follows  the  general  rule  that  "like  dis¬ 
solves  like".  Hence  we  may  also  say  that 

5.  "Like  permeates  through  like." 

6.  Since  polyethylene  is  itself  relatively  non-polar,  P-Factors  in¬ 
variably  decrease  as  the  polarity  of  the  penetrant  increases. 

7.  There  are  several  reliable  ways  of  indicating  the  polarity  of  a 
penetrant. 

8.  Qualitatively  at  least 

a.  Dunkel's  group  contributions  to  AH^p  and 

b.  the  Cohesive  Energy  Density  of  the  penetrant  are 
fairly  reliable  guides  to  polarity. 

9.  Quantitatively  it  has  been  found  that  the  so-called 

a.  Polarity  Index  AAH  as  determined  from  experimental 
heats  of  vaporization  and  the  so-called 

b.  Flory-Huggins  Polymer  Interaction  Parameter 
as  determined  by  swelling  equilibrium 

are  better  guides  to  polarity  in  polyethylene  -  organic  solvent  systems. 

10.  A  correlation  was  found  between  the  Polarity  Index  AAH  and  the 
temperature  variation  of  the  Flory-Huggins  parameter. 

11.  An  oxygen  atom  may  often  be  a  seat  of  polarity  in  a  molecule.  It 
is,  however,  very  important  to  distinguish  the  type  of  linkage. 

12.  All  of  the  indices  mentioned  in  8  and  9  agree  in  assigning  the 
greatest  polarity  to  the  -OH  or  hydroxyl  group,  which  can  hydrogen  bond  with 
itself,  the  next  greatest  polarity  for  oxygen  linkages  to  the  C  =  O  or  carbonyl 
group,  and  the  least  polarity  to  the  C-O-C  linkage.  Thus 

13.  The  alcohols  were  the  most  polar  and  the  least  permeable  of  all 
the  substances  studied.  The  esters  and  ketones  were  always  somewhat  more 
permeable,  and  the  ethers  resembled  the  corresponding  hydrocarbons. 

14.  The  carboxylic  acids  are  more  permeable  than  their  polarity 
indexes  would  suggest.  They  are,  in  fact,  effectively  less  polar  than  the  cor¬ 
responding  alcohols. 

15.  Considerable  evidence  points  to  the  existence  of  carboxylic  acid 
dimers  in  polethylene.  This  largely  explains  their  seemingly  anomolous 
behavior. 
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16.  Dimerization  tendencies  also  explain  why  the  heat  of  solution  of 
penetrant  in  polymer  is  generally  significantly  less  thanAAH. 

17.  In  general,  though,  any  hydrogen  bonded  substance  will  have  a  very 
low  permeability  factor. 

18.  With  regard  also  to  the  polarity  factor,  there  is  evidence  for  in¬ 
ternal  shielding  effects  in  the  higher  aliphatic  compounds,  which  operates  to 
increase  S  and  P  slightly. 

19.  This  internal  shielding  effect  also  explains  why  the  order  of  S  is 
the  reverse  of  Ep  for  members  within  a  single  homologous  series,  except  in 
the  extreme  cases  of  the  hydrocarbons  and  the  alcohols,  which  can  be  ex¬ 
plained. 

20.  The  limited  swelling  that  polyethylene  undergoes  in  even  the  best 
solvents  at  room  temperature  resides  in  its  partial  crystallinity.  Otherwise 
it  would  dissolve  in  most  hydrocarbons. 

21.  The  partial  crystallinity  of  polyethylene,  however,  gives  rise  to 
an  "anti-freeze  effect"  which  seriously  complicates  the  quantitative  develop¬ 
ment  of  the  theory  of  permeability. 

22.  It  has  been  assumed  that  the  penetrants  do  not  penetrate  the 
crystalline  regions  in  polyethylene. 

23.  The  solubility  coefficient,  however,  is  not  the  whole  story. 

a.  Although  S  always  increases  with  temperature,  the 
P-Factor  itself  increases  much  more  rapidly. 

b.  Indeed,  Ep  was  in  every  instance  greater  than  12 
whereas  Eg  was  always  less  than  5. 

c.  Furthermore,  S  is  not  very  much  dependent  upon 
molecular  geometry,  but  P  very  definitely  is. 

d.  Finally,  the  effect  of  polarity  extends  beyond  S. 

24.  A  quantitative  picture  of  the  permeability  of  polyethylene  to  organic 
solvents  is  best  based  upon  the  relation: 

p  -  p  e-E_/RT. 
o  ^ 

25.  The  Master  logP  vs.  1 /T  Plot  can  be  used  to 

a.  Determine  the  parameters  E  and  Pq;  to 

b.  Check  the  over-all  consistency  of  the  data;  to 

c.  Smooth  the  experimental  data;  to 

d.  Extrapolate  the  experimental  data;  and  to 

e.  Interpolate  with  the  experimental  data  between 

i.  two  temperatures,  same  compound,  or 
ii.  two  compounds,  same  temperature. 


WADC  TR  53-133  Pt  2 


65 


26.  The  size  and  shape  of  the  penetrant  affect  P  primarily  through 
the  kinetic  factor  D. 

27.  It  has  been  proposed  to  use  the  molar  volume  as  an  index  of  size 
and  the  molar  volume  divided  by  the  length  of  the  corresponding  Fisher- 
Hirschfelder- Taylor  molecular  model  as  an  index  of  shape. 

28.  The  diffusion  coefficient  is  also  significantly  dependent  upon  the 
polarity  index  AAH. 

29.  Quantitatively,  we  can  say  that 

a.  E  increases  0.0348  Kcal/mole  per  cc  of  molar  volume; 

b.  Ep  increases  0.75  Kcal/mole  per  unit  increase 
in  the  shape  factor  V/L;  and 

c.  Ep  increases  2.4  Kcal/mole  per  Kcal  of  polarity  as 
measured  by  AAH;  or  altogether, 

d.  Ep  =  0.0348V-#-  0.  75V /L  2.4A&H. 

30.  In  some  cases  it  is  necessary  to  correct  29d  as  follows: 

a.  Subtract  1  for  para- substituted  compounds. 

b.  Subtract  2  for  all  alcohols. 

c.  Add  1.3  for  all  aldehydes  {This  rule  is  based  on 
data  for  only  two  compounds). 

31.  At  the  moment  AAH  for  the  carboxylic  acids  must  be  determined 
from  swelling  measurements. 

32.  In  general  E  is  more  important  than  Pq  in  determining  relative 
P  values  at  room  temperature. 

33.  A  rather  remarkable  correlation  between  E  and  logPQ  has  been 
found  and  partially  explained. 

34.  Using  this  logP  vs.  E  Plot, 

o  p 

a.  A  single  P- Factor  determination,  say  at  room 
temperature,  suffices  to  determine  the  permeability 
behavior  at  all  other  temperatures; 

also 

b.  an  absolute  value  for  the  P- Factor  can  be 
calculated  at  any  temperature  given  only  E^,  which 
we  can  calculate. 


That  is  to  say,  it  is  only  necessary  to  have  the  usual  vital  statistics 
on  the  compound  in  question,  namely  its  density,  molecular  weight,  model 
length,  and  heat  of  vaporization,  and  one  can  calculate  good  values  for  its 
permeability  factor. 
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The  value  of  the  equation; 


P  =  P  e"Ep/RT 
o  r 

may  be  observed  in  Table  2 1  in  which  the  P-factors  calculated  by  the  above 
equation  are  compared  with  observed  values.  In  Table  21  are  listed  most  of 
the  substances  examined  in  this  report  (actually  all  those  given  in  Table  13). 
The  agreement  between  the  fourth  and  fifth  columns  is  fairly  impressive.  As 
indicated  by  the  last  column,  the  calculated  values  are  clearly  of  the  right 
order  of  magnitude  (off  by  less  than  a  factor  of  10),  and,  with  one  exception, 
better  even  than  the  hoped  for  factor  of  2.  Indeed,  in  50%  of  the  cases  the 
calculated  and  observed  values  agree  to  within  the  limit  of  the  probable  ex¬ 
perimental  error. 
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Table  21 


A  FINAL  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT: 


The  Prediction  of  Permeability  Factors  at  70°F  for 
the  Diffusion  of  Organic  Solvents  through  Polyethylene. 


Substance 

caka 

P 

b 

loSPo 

_calcc 

y0°F 

v  d 

pobs 

70  °F 

pcalc  pob 
70°F7  70 

n-Octanol 

23.  5 

17.  17 

0.  52 

0.  50 

1. 05 

ter -Butanol 

26.  7 

19.  36 

0.  34 

0.26 

1.  30 

sec -Butanol 

22.  7 

16.  62 

0.  54 

0.  60 

0.  88 

n-Butanol 

22.4 

16.40 

0.  57 

0.46 

1.25 

n-Propanol 

22.4 

16.40 

0.  57 

0.40 

1.44 

Phenol e 

20.  9 

15.  38 

0.  72 

0.48 

1.44 

Butyr  aldehyde  e 

19.  0 

15.  10 

9.  77 

10.  13 

0.97 

Benzaldehyde  e 

19.2 

15.20 

8.  71 

6.  80 

1.28 

Acetone  e 

18.  3 

14.76 

14.8 

6.  75 

2.  19 

Methyl  Ethyl  Ketone  6 

18.  3 

14.76 

14.  8 

12.  6 

1.  16 

Ethyl  Acetate  e 

18.  1 

14.  66 

16.  6 

17.0 

0.  98 

Amyl  Acetate 

19.  5 

15.  35 

7.41 

8.  30 

0.  89 

Dibutyl  Ether  e 

16.  3 

13.  80 

49 

88 

0.  55 

Pentane 

13.0 

12.22 

371 

526 

0.  70 

Heptane 

14.0 

12.  68 

190 

270 

0.  70 

Decane 

15.  7 

13.  50 

69.  2 

71.2 

0.  97 

Tetradecane 

18.  0 

14.  60 

17.  0 

17.0 

1. 00 

Benzene 

13.  3 

12.  52 

436 

440 

0.99 

o- Xylene 

16.  0 

14.  38 

316 

290 

1.  09 

p- Xylene 

13.  6 

12.  73 

427 

500 

0.  85 

Chlorobenzene 

14.  3 

13.22 

398 

500 

0.  79 

p  -  Chlor  ot  olue  ne 

16.2 

14.  51 

302 

300 

1. 01 

6 

Cyclohexane 

16.  8 

14.  93 

282 

285 

0.99 

Dipentene  * 

15.9 

13.  95 

138 

128 

1. 08 

Ep  expressed  as  kcal/mol 
P  r  grams/24  hrs/0.001  in,  / 100  in.  ^ 
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Notes  for  Table  21 


a  Taken  from  Table  13. 

b  Obtained  using  the  values  in  the  second  column  and  Figure 

c  Calculated  using  equations  (l-2a)  and  (1-3). 

d  Smoothed  experimental  values  obtained  from  the  Master  Plot, 

e  See  section  25. 

f  See  section  26,  point  (1). 
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A  few  sample  calculations  of  the  P-factor  using  the  Bent  Equation  and  ap¬ 
pearing  in  Table  21  are  included. 

n-Octanol 

P  =  P  e_EP/RT  or 
o 

log  P  =  -Ep/RT  4-  logPQ 

=  -  23.5/  (2.  3x  1. 987x10"3x  294)  +  17.17 
=  -  23.  5/  (1.  345)  -t-  17.  17 
=  -  17.45  +  17.  17 
=  -0.28  or  9.72-10 

P  =0.52  gms/24hrs/0.  OOlin/lOOin^ 

Tetradecane 

log  P  =  -Ep/RT  -+-  log  PQ 

=  -18.0/  (2.3  X  1 . 987xl0"3x  294)  +  14.60 
=  _  18.  0/  (1.  345)  +  14.  60 
=  -13.  37  -f  14.60 
=  1.23 

P  z  17.0  gms /24hr  s /0.  001in/100  in^ 

Benzene 

log  P  =  -Ep/RT  +  log  PQ 

=  -13.3/  {2.3  k  1.987*103X  294)4-12.52 
i  -13.  3/(1.  345)  +  12.  52 
=  -9.  88  4-  12.  52 
=  2.64 

P  =  436  gms /24hr  s /0.  001  in/100  in^ 
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MASTER  PLOT  P  «S  y 

'  & 


100,000. 


10,000. 


1,000. 


100. 


10. 


p  XYLENE _ 

t«r  BUTANOL  _ 


BENZENE 


AMYL  ACETATE 


DECANE 


TCTRADECANE 


METHYL  ETHYL  KETONE 


D  OCTYL  ALCOHOL _ 


n  BUTYRIC  ACID 


n  PROPYL  ALCOHOL 


ACETONE 


ACETIC  ANHYDRIDE 


BENZOIC  ACID 


DIACETONE  ALCOHOL* 


NITROETHANE 


ACETIC  ACID  * 


FORMIC  ACID 


WATER 


- 

1 

10. 

1.0 

x 

n-PENTANE 

p-XYLENE  jf 

BENZENE  X  X^ 


rt-HEPTANE  _ _ 


CHLOROBENZENE  OX  0.| 


rs7 


p-CHLOROTOLUENE  -+-. 

o-XYLENE _ / 

CYCLOHEXANE  _ 


decane _ ®/o.oi 


ACET^amc  '^VX>00t  4.3  4.2  4.1  4.0  3.9  3.8  3.7  3.6  3.5  3.4  3.3 

tetrapecane  o// 

ETHYL  ACETATE  /  -4  0.  -3  1.  -2  0.  -9.4  +1.4  +14.  +2  7.  +4  1.  +5  9.  +7  0.  +8  6. 


3.4  3.3  Z.Z  3.1  3.0  2.9  2.8  2.7  2.6  2.5  2.4 

+  70.  +86.  +I04.  +  I22.  +  I4I.  +162.  +  r83,  +  207.+234.  +  26l.  +289. 


This  is  a  'selected'  bibliography.  It  is  not  intended  to  be  complete. 
It  should  nevertheless  afford  an  easy  entry  to  most  of  the  pertinent  litera¬ 
ture  on  permeability,  solubility,  and  diffusion  in  polymer  systems. 


I.  ARTICLES. 


A.  Permeability  in  General. 


Van  Amerongen  appears  to  be  nearly  the  only  investigator  who  has 
published  work  concerned  primarily  with  the  permeability  factor  itself.  Two 
of  his  articles  are  listed  below. 

1.  "The  Permeability  of  Different  Rubbers  to  Gases  and  Its  Relation 
to  Diffusivity  and  Solubility.  " 

G.  J.  Van  Amerongen,  J^.  App.  Phys.  17,  972  (1946). 

A  very  pertinent  article.  Many  of  his  conclusions  are 
similar  to  those  reached  in  this  report. 

2.  "Solubility,  Diffusion,  and  Permeation  of  Gases  in  Guttapercha.  " 

G.  J.  Van  Amerongen,  Poly.  Sci.  2,  381  (1947). 

A  typically  readable  Van  Amerongen  article. 

B .  Diffusion. 

Although  the  literature  on  this  subject  does  not  include  much  on  poly¬ 
ethylene,  diffusion  in  polymers  has  received  considerable  attention,  particu¬ 
larly  by  a  group  of  investigators  at  ICI.  Their  work  is  generally  published 
in  the  Transactions  of  the  F  ar  ad  ay  Society. 

1.  "The  Diffusion  of  Some  Organic  Substances  in  Polystyrene.  " 

G.  S.  Park,  Trans.  Far .  Soc .  47,  1007  (1951). 

An  interesting  paper.  Attempts  to  describe  quantitatively 
the  effect  of  size  and  shape  on  the  diffusion  constant  (The  'Smallest- 
Hole"  Theory). 

2.  "The  Diffusion  of  Some  Halo-Methanes  in  Polystyrene." 

G.  S.  Park,  Trans.  Far.  Soc.,  46,  684  (1950). 
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B.  Diffusion  cont. 


The  predecessor  to  the  paper  above.  Considers  only 
the  molar  volume  of  the  penetrant.  Also  gives  a  brief  discussion 
of  the  diffusion  coefficient  in  terms  of  the  transition  state  theory. 

3.  "Diffusion  in  High  Polymers:  Some  Anomalies  and  Their 
Significance.  " 

J.  Crank  and  G.  S.  Park,  Trans.  Far.  Soc.  51,  10072  (1951). 

Anomalies  are  attributed  to  (a)  variable  surface  concen¬ 
tration,  (b)  diffusion  coefficient  dependence  on  the  history  of  the 
diffusion  process,  and  (c)  stresses  exerted  on  one  part  of  the 
polymer  sheet  by  other  parts. 

4.  "Diffusion  of  Organic  Vapors  into  Polyvinyl  Acetate.  " 

R.  J.  Kokes  and  F.  A.  Long,  J.  Am.  Chem.  Soc .  75,  6142  (1953). 

A  report  of  some  interesting  recent  work  with  many 
observations  that  are  of  significance  for  polyethylene.  The  influ¬ 
ence  of  penetrant  concentration,  size,  and  polarity  upon  D  is 
noted. 

5.  "Theory  of  DQ  for  Atomic  Diffusion  in  Metals." 

Clarence  Zener,  _L  App.  Phys.  22,  372  (1951). 

Included  here  for  its  clarifying  discussion  of  the  relation 

of  S*  to  D  . 

o 

C.  Solubility. 

The  theoretical  work  in  the  high  polymer  field  of  solubility  was  iniati- 
ated  essentially  independently  in  1942  by  Flory  and  Huggins.  E.  A. 
Guggenheim  has  recently  thoroughly  reviewed  the  theory  in  a  book,  "Mixtures". 
These  treatments  are  all  quite  mathematical.  The  emphasis  in  the  references 
listed  below  is  by  contrast  largely  experimental. 

1.  "Influence  of  Structure  on  Polymer -Liquid  Interaction. 

I.  Relative  and  Absolute  Values  of  Swelling  Equilibria.  " 

G.  Salomon  and  G.  J.  Van  Amerongen,  £.  Poly.  Sci.  2_,  355  (1947). 

A  good  introductory  discussion;  possibly  the  most 
readable  in  the  literature  on  this  somewhat  difficult  subject. 
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C.  Solubility  cont. 


2.  "Influence  of  Structure  on  Polymer -Liquid  Interaction. 

III.  Swelling  and  Mechanical  Properties  of  Some  Partly 
Crystallized  Polymers." 

G.  Salomon,  J.  Poly.  Sci.  3,  776  (1948). 

Continues  along  the  lines  of  the  previous  reference. 

3.  "The  Interaction  Between  Rubber  and  Liquids.  VIII.  A  New 
Examination  of  the  Thermodynamic  Properties  of  the  System  Rubber 

t Benzene.  " 

Geoffrey  Gee  and  W.  J.  C.  Orr,  Trans.  Far.  Soc.  42,  507  (1946). 

Suggests  .a  number  of  experiments  that  might  be  tried  on 
polyethylene . 

4.  "The  Phase  Equilibria  Between  A  Crystalline  Polymer  and  Solvents. 

I.  The  Effect  of  Polymer  Chain  Length  on  the  Solubility  and 
Swelling  of  Polythene. 

II.  The  Effect  of  Solvent  Type  on  the  Solubility  and  Swelling  of 
Polythene.  " 

R.  B.  Richards,  Trans.  F ar .  Soc.  42,  10-28  (1946). 

Highly  recommended  reading  for  anyone  interested  in 
polyethylene.  These  articles  contain  a  lot  of  information  that 
might  profitably  be  correlated  with  that  of  the  present  study. 

5.  "Swelling  and  Solubility  in  Mixed  Liquids." 

Geoffrey  Gee,  Trans ,  Far .  Soc.  40,  468  (1944). 

More  good  reading.  Use  is  made  of  the  cohesive  energy 
density  argument. 

6.  "Determination  of  Polymer -Liquid  Interaction  by  Swelling  measure¬ 
ments.  " 

Paul  Doty  and  Helen  S.  Zable,  J.  Poly.  Sci.  1 ,  90  (1946). 

A  fairly  detailed  report  on  the  swelling  of  a  cross- 
linked  polymer  in  a  wide  variety  of  liquids.  Includes  a  relatively 
non-mathematical  discussion  of  Flory-Huggins  interaction  parameter. 

7.  "Some  Thermodynamic  Properties  of  Slightly  Cross-Linked  Styrene  - 
Divinylbenzene  Gels." 

R.  F.  Boyer  and  R.  S.  Spencer,  J.  Poly.  Sci.  3,  97  (1948). 
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C.  Solubility  cont. 


Very  similar  in  content  to  the  previous  reference. 
Quite  complete  so  far  as  the  experimental  work  goes  with  many 
references. 

8.  "Thermodynamics  of  High-Polymer  Solutions.  III.  Swelling  of 
Cross-Linked  Rubber." 

R.  L.  Scott  and  M.  Magat,  _L  Poly.  Sci.  555  (1949). 

A  description  of  some  absorption  experiments  that  are 
very  similar  to  those  of  the  present  project. 

9.  "Equilibrium  Sorption  of  Several  Organic  Diluents  in  Polyvinyl 
Acetate.  " 

R.  J.  Kokes,  A.  R.  DiPietro  and  F.  A.  Long,  J.  Am.  Chem. 
Soc.  75,  6319  (1953).  "  ~ 

Another  recent  article  by  F.  A.  Long  and  co-workers. 
Includes  an  application  of  equation  (18-9). 


D.  The  Physical  Properties  of  Polyethylene. 

One  of  the  distinguishing  features  of  polyethylene  is  that  it  is  a  partial¬ 
ly  crystalline  at  room  temperature.  The  evidence  for  this  consists  primarily 
of  the  pioneer  work  of  the  ICI  research  team  on  the  heat  capacity,  heat  of 
solution,  density,  and  x-ray  diffraction  of  polyethylene.  More  recently  a 
group  at  Du  Pont  has  made  some  contributions  to  this  field. 

1.  "The  Effect  of  Temperature  on  the  Density  of  Polyethylene." 

E.  Hunter  and  W.  G.  Oakes,  Trans .  Far.  Soc.  41,  49  (1945). 

Perhaps  the  simplest  way  of  following  the  change  in 
crystallinity  with  temperature. 

2.  "The  Heat  Capacity,  Heat  of  Solution,  and  Crystallinity  of  Poly¬ 
ethylene  .  " 

H.  C.  Raine,  R.  B.  Richards  and  H.  Ryder,  Trans.  Far.  Soc. 

41,  56  (1945). 

Contains  data  relating  to  the  "anti-freeze”  effect. 

3.  "The  Melting  of  Polythene." 

R.  B.  Richards,  Trans.  Far.  Soc.  41,  127  (1945). 

A  noteworthy  attempt  to  give  a  theoretical  explanation  of 
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D.  The  Physical  Properties  of  Polyethylene,  cont. 


the  phase  relationships  in  pure  polyethylene  and  swollen  poly¬ 
ethylene;  is  successful  in  only  a  semi-quantitative  way.  The  dis¬ 
cussion  is  sometimes  a  bit  mathematical.  (Supersedes  Frith  and 
Tuckett's  treatment,  "The  Melting  of  Crystalline  Polymers", 

Trans.  Far.  Soc.  40,  251  (1944)  ) 

4.  "The  Texture  of  Polythene.  " 

C.  W.  Bunn  and  T.  C.  Alcock,  Trans.  Far.  Soc.  41 ,  317  (1945). 

Beautiful  evidence  for  the  crystallinity  of  polyethylene. 

5.  "The  Molecular  Structure  of  Polyethylene.  I.  Chain  Branching  in 
Polyethylene  During  Polymerization.  " 

M.  J.  Roedel,  JL  Am.  Chem.  Soc .  75,  6110  (1953). 

This  is  the  first  of  a  series  of  five  articles  by  the 
Du  Pont  group  that  appear  consecutively  in  the  same  issue  of 
JACS.  We  mention  them  here  primarily  because  they  are  of  some 
interest  and  afford  a  very  recent  entry  to  the  literature  in  this 
field.  The  titles  of  the  remaining  articles  are  given  below. 

II.  Determination  of  Short  Chain  Branching. 

III.  Determination  of  Long  Chain  Branching. 

IV.  Kinetic  Calculations  of  the  Effect  of  Branching  on 
Molecular  Weight  Distribution. 

V.  The  Effect  of  Chain  Branching  and  Molecular  Weight 
on  Physical  Properties. 

6.  "The  Mechanical  Properties  of  High  Polymers.  " 

H.  Mark,  Trans.  Far.  Soc.  43,  447  (1947). 

A  lecture  by  Professor  Mark  which  nicely  surveys  the 
polymer  field. 


n.  BOOKS. 

There  are  several  books  that  have  been  found  useful  in  this  problem. 
Listed  in  order  of  more  or  less  decreasing  interest,  they  are: 

1.  "Elastomers  and  Plastomers.  I.  General  Theory." 

Edited  by  R.  Houwink.  Elsevier  Publishing  Co.  (1950) 

The  chapters  "Molecular  Constitution",  "Physics  and 
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II.  BOOKS  cont. 


Structure"  (Section  5  is  on  Permeability),  and  "Polymer -Liquid 
Interaction"  by  Van  Amerongen  form  an  excellent  introduction  to 
the  polyethylene  permeability  problem. 

2.  "The  Solubility  of  Nonelectrolytes.  " 

J.  H.  Hildebrand  and  R.  L.  Scott.  Reinhold.  (1950) 

The  new  third  edition  of  this  classic  contains  a  chapter 
on  high  polmer  solutions  with  seventy-one  references  to  the 
literature. 

3.  "Advances  in  Colloid  Science."  Vol.  II. 

Edited  by  H.  Mark  and  G.  S.  Whitby.  Inter  science .  (1946) 

Gee  contributed  a  readable  chapter  to  this  volume  on 
"The  Thermodynamic  Study  of  Rubber  Solutions  and  Gels." 

4.  "Diffusion  in  and  through  Solids." 

R.  M.  Barrer.  Cambridge  University  Press.  (1951) 

This  book  which  first  appeared  in  1941  contains  much 
valuable  data  on  diffusion. 

Finally,  mention  should  be  made  of  the  invaluable  "Annual  Reviews  of 
Physical  Chemistry"  which  was  established  in  1950.  This  review  almost 
always  contains  something  on  polymers.  Volume  4,  for  example,  has  a 
chapter  on  the  "Physical  Properties  of  High  Polymers"  with  310  citations  to 
the  literature,  and  volume  5  is  also  scheduled  to  have  a  chapter  on  polymers 
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APPENDIX  A 


Derivation  of  the  Ideal  Solubility  Equation:  N  -  expl-AH^^/RT) 


'  According  to  classical  thermodynamics,  when  a  polymer  is  saturated 
with  a  solvent  designated  here  as  component  2, 

m-  p  sy^L 
2  2 

p 

where^fjL-2  =  the  so-called  chemical  potential  or  the  partial  molal  free  energy 
of  the  solvent  in  the  saturated  polymer  andyu.^  =  the  chemical  potential  of  the 
solvent  in  the  pure  solvent  phase.  Furthermore,  we  have  that 


■  ?  =  H?  -  TS?  and 


where 


jul  L  s  HL  -  TSL 
^  2  2  2 

AH  solution  a  -  H^1, 

SL  «  kin  WL,  and 

2  2 


S?  *  3  S P  ;  SP  »  kin  WP 

3-2 


p 

at  this  point  it  becomes  necessary  to  make  some  assumptions  regarding  W 
We  shall  assume  that  WP  can  be  written  as 


wp  - 

Wp  WP  Wc 

1  2 

where 

II 

(WL  j  n2, 

Wc  - 

(— sites )  * 

(-sites  ”-2^  •  n2  • 

and  where  WP  is  independent  of  n ^  as  is  ngjtes  the  total  number  of  positions 
in  the  polymer  available  to  organic  solvent  molecules.  The  equations  above 
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essentially  define  what  we  have  called  an  ideal  solution. 


Now,  since  S  ^  =  klnW  ^ 

=  k  JinW^  +  lnW^t  lnW^j 

and  since 

Sin  W? 

w 

0n2 

'  lnW2  • 

^^soln  = 

T(S^  -  sh  *  T{klnW^  -klnwj"  -k  d  lnW( 

n2 

= 

RT  ln  nsites“n2  ~,RT  In 

nsites 

n2 

n2 

for  dilute  solutions  where  n^  < 

<  n  . , 

^  sites 

calling  n2/nsites  =  N*  *kis  6ives  N  =  e  “J^Lsoln 
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APPENDIX  B 


The  Heats  of  Solution  of  Liquids  in  Polyethylene. 

Let  us  return  to  our  earlier  example  and  consider  the  mechanism  by 
which  a  molecule  of  benzene  dissolves  in  polyethylene. 

First  it  is  necessary  to  separate  it  from  the  pure  liquid  benzene  phase 
by  breaking  the  intermolecular  bonds  that  it  forms  with  neighboring  benzene 
molecules.  If  we  were  to  do  this  for  a  whole  mole  of  molecules  we  would 
write 

(17.  1)  Benzene  (liquid)  *  Benzene(gas), 

and  the  heat  of  this  reaction  is  just  the  heat  of  vaporization  of  pure  benzene. 
It  is  symbolized  by  &Hva_orization  which  stands  for  the  difference  Hbenzene 

as  a  gas  “  ^benzene  as  a  liquid' 

The  second  step  in  the  solution  process  is  then  to  place  this  free 
(gaseous)  benzene  molecule  in  the  polymer, 

(17.2)  Benzene  (gas)  *  Benzene  (polymer). 

The  sum  of  reactions  (17.1)  and  (17.2)  is  reaction  (13)  which  may  be  ab¬ 
breviated  as 

(17.3)  Benzene  (liquid)  =  Benz  ene  (polymer) . 

Calling  the  heat  of  the  reaction  (17.2)  ^H§g^ution,  we  have  that  the  heat  of 
the  net  reaction,  equation  (17.3),  is  just  the  sum  of  the  heats  of  the  two 
reactions  (17.1)  and  (17.2).  That  is, 

{*8)  Absolution  =  ^bv.ap0rization+  AH®oju^on. 

Thus  we  axe  led  from  the  permeability  P  to  the  solubility  S  (P  :  DS) 
to  the  heat  of  solution  AHgolution  (N  =<" AHsoln/RT»  and  finally  to  the  heat 

of  vaporization  and  AH® (equation  (18)  ). 

Of  these  latter  two  terms,  AH^gp  is  of  course  always  positive  whereas 
Solution  has  just  the  opposite  sign  principally  because  of  the  intermolecular 
attractive  forces  which  play  a  relatively  important  role  in  condensed  phases 
but  are  unable  to  exert  their  full  effect  in  the  comparatively  dilute  gaseous 
phase.  Since  it  is  to  be  expected  that  AHsojn  will  be  positive,  the  absolute 
magnitude  of  must  be  less  than  that  of  AHyap 
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The  interplay  of  these  three  quantities  can  perhaps  best  be  discussed 
by  considering  two  contrasting  cases:  propanol  and  decane. 

It  has  been  established  that  both  the  permeability  and  the  solubility  {in 
polyethylene)  of  decane, 

CH3CH2CH2CH2CH2CH2CH2CH2CH2CH3,  are  significantly  greater  than  those 
of  propanol,  CH3CH2CH2OH.  From  equation .  (15),  N  8  (“^^soln^*^,  we 
conclude  therefore  that  the  heat  of  solution  of  decane  in  polyethylene  is  less 
than  the  heat  of  solution  of  propanol  in  polyethylene.  Looking  now  at  equa¬ 
tion  (18),  we  might  be  inclined  initially  to  conclude  that  the  heat  of  vapor¬ 
ization  of  decane  is  significantly  less  than  the  heat  of  vaporization  of 
propanol,  yet,  in  fact,  they  are  nearly  equal  (we  have  already  seen  that  the 
heat  of  vaporization  of  propanol  is  about  11.6  kcal/mole.  Using  table  4  it 
is  estimated  that  the  heat  of  vaporization  of  decane  should  be  8(1.0)  +■  2(1.8) 

*  11.6  kcal/mole.). 

The  reason  for  this  near  equality  and  actually  relatively  high  heat  of 
vaporization  is  that  in  the  case  of  decane  although  the  intermolecular  force 
field  surrounding  the  molecule  is  at  no  point  very  intense  it  does  operate 
over  a  large  area  and  manages  to  accumulate  to  11.6  kcal/mole;  and  in  the 
case  of  the  much  smaller  propanol  molecule  hydrogen  bonds  are  largely 
responsible  for  this  high  heat  of  vaporization.  Indeed,  from  Table  4  one 
-OH  group  is  in  this  respect  equivalent  to  roughly  7  -CH2-  groups. 

Now  when  we  come  to  consider  the  situation  to  be  expected  in  the 
polymer,  we  see  that  a  decane  molecule  experiences  very  much  the  same 
general  dispersion  forces  that  exist  in  pure  decane,  which  implies  that 
^^lofn  i*1  this  case  be  comparable  to  AHygp.  However,  for  propanol 

the  situation  is  quite  different  since  polyethylene  is  a  nonpolar,  essentially 
neutral  or  non-specific  solvent  medium.  Thus  the  0H*^0  bonds  that  play  such 
an  important  role  in  pure  propanol  are  largely  non-existent  in  the  poly¬ 
ethylene  system,  causing  AHSas  to  subtract  off  much  less  from  AHyaD  than 
in  the  case  of  decane  (equation  (18)  ).  This  then  explains  why  the  heat  of 
solution  of  propanol  in  polyethylene  can  be  expected  to  be  significantly 
greater  than  that  of  decane  even  though  both  have  essentially  the  same  heat 
of  vaporization. 

More  generally,  since  the  -CH2-  groups  in  polyethylene  have  a  rela¬ 
tively  small  polarizability  (see  C-H  in  Table  V)  and  dipole  moment  (see 
HC-  in  Table  III)  so  that  London-type  and  dipole -dipole-  type  interactions 
with  groups  on  other  molecules  will  at  most  be  quite  non-specific  and  small 
and  furthermore  since  -CH2-  groups  do  not  form  hydrogen  bonds,  it  is  to  be 
expected  that  whenever  these  specific  intermolecular  interactions  do  exist  in 
the  pure  solvent  phase  (as  for  so-called  polar  liquids),  ^Hvap0rjzati0n  will 
be  large  compared  to  AH|a|n>  and,  therefore  ^^so|utjon  will  in  suchinstances 
be  large  when  AHvaporization  itself  is  large. 
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